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ABSTRACT 
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The chemistry of gold dissolution in alkaline cyanide solution has continually 
received attention and new rate equations expressing the gold leaching are still 
developed. The effect of leaching parameters on gold cyanidation is studied in this 
work in order to optimize the leaching process. A gold leaching model, based on 
the well-known shrinking-core model, is presented in this work. It is proposed that 
the reaction takes place at the reacting particle surface which is continuously 
reduced as the reaction proceeds. The model parameters are estimated by 
comparing experimental data and simulations. The experimental data used in this 
work was obtained from Ling et al. (1996) and de Andrade Lima and Hodouin 
(2005). Two different rate equations, where the unreacted amount of gold is 
considered in one equation, are investigated. In this work, it is presented that the 
reaction at the surface is the rate controlling step since there is no internal diffusion 
limitation. The model considering the effect of non-reacting gold shows that the 
reaction orders are consistent with the experimental observations reported by Ling 
et al. (1996) and de Andrade Lima and Hodouin (2005). However, it should be 
noted that the model obtained in this work is based on assumptions of no side 
reactions, no solid-liquid mass transfer resistances and no effect from temperature. 
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NOMENCLATURE 
Roman Letters 



A surface area of gold [m 2 ] 

A i surface area of one particle during leaching [m 2 ] 

A r reacting surface area of particles [m 2 ] 

A n inner surface area of passive layer [m 2 ] 

A r outer surface area of passive layer [m ] 

a specific interfacial area [m 2 m R 3 ] 

b Freundlich adsorption constant [-] 

C concentration in spherical diffusion layer [mol m L " 3 ] 

C Au gold concentration in particles [mol m R " 3 ] 

C Au ini gold initial concentration in particles [mol m R 3 ] 

C°° Au gold concentration in particles after infinite leaching time [mol m R 3 ] 

"3 

Cau,c S°W concentration in carbon [mol m" ] 

C'Z,l concentration in solution at initial time [mol m L " ] 

Cau,l g 0 ^ concentration in solution [mol m L " 3 ] 

C Au L (t) gold concentration in solution at time t [mol m L " 3 ] 

Ca u ,s surface concentration of gold in particles [mol m" ] 

C CA ,_ cyanide concentration [mol m L " 3 ] 

C, concentration in solution at inner passive layer [mol m L " 3 ] 

C 0 concentration in solution at outer passive layer [mol m L " 3 ] 
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"3 

C Qi oxygen concentration [mol m L " ] 

C 0i S dissolved oxygen concentration at the reacting surface [mol m L " ] 
Cj , c 2 coefficients in Equation (27) 
c 3 , c 4 coefficients in Equation (30) 

2 1 

D CN . diffusion coefficient of cyanide [m s" ] 

D e effective diffusivity in particles [m 2 s" 1 ] 

2 1 

D m molecular diffusivity [m s" ] 

2 1 

D 0i diffusion coefficient of oxygen [m s" ] 

d average diameter of particles [m] 

Jd impeller diameter [m] 

dp particle diameter [m] 

Jt tank diameter [m] 

E° standard reaction potential [V] 

E h reduction potential [V] 

F the Faraday constant (96,485.34) [C mol" 1 ] 

Ga Galileo number, dp 3 Apgpi 2 //j [-] 

g gravitational acceleration [m s" 2 ] 

H distance traveled by eddy in concentrated boundary 

layer surrounding the particle [m] 

k overall rate constant in Equation (25) and (26) 



k l overall rate constant in Equations (42), (52) and (54) [(mol/m 3 )" 2 ' 35 h" , 

(mol/m 3 ) 0 - 668 !!- 1 ] 

k 2 overall rate constant in Equations (43), (53) and (55) [(mol/m 3 )" 2 ' 731 h" , 

(mol/m 3 )' 2 - 011 h" 1 ] 
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k a rate constant for anodic half-reaction [C m mol" 1 s" 1 ] 

k ads rate constant for gold adsorption [s" 1 ] 

k CN - rate constant in Equation (23) [s" 1 ] 

k c rate, constant for cathodic half-reaction [C m mol" 1 s" 1 ] 

k L liquid-gas mass transfer coefficient [m s" 1 ] 

k L a volumetric mass transfer coefficient [s" 1 ] 

k 0i rate constant in Equation (24) [s" 1 ] 

&sl solid-liquid mass transfer coefficient [m s" 1 ] 

m constant in Equation (33) [-] 

m T total mass of particles [kg] 

./V stirrer speed [s" 1 ] 

N 02 diffusion rate of oxygen [mol s" 1 ] 

Nau amount of gold in particles [mol] 

iVpo power number, P/pN 3 dn 5 [-] 

n number of particles per volume [m R " ] 

n' Freundlich exponent [-] 

P power input under gassed conditions [W] 

Re^ Reynolds number based on impeller diameter, Ndvplp [-] 

Re-i Reynolds number based on tank diameter, NT 2 pin [-] 

Rep Reynolds number based on particle diameter, dpdjconln [-] 

Re E Reynolds number based on Kolmogoroffs theory of 

Isotropic turbulence, (sdp 4 /v 3 ) lB [-] 

r radial coordinate [m] 

3 1 

r CN _ reaction rate for cyanide consumption in Equation (29) [kg m" s" ] 

r, inner radius of passive layer [m] 
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r„ outer radius of passive layer [m] 

r s reaction rate at surface [mol m" 2 s" 1 ] 

Sc Schmidt number, plpD m [-] 

Shp Sherwood number based on particle diameter, k S hdp/D m [-] 

Shi Sherwood number based on tank diameter, ksLdjlD m [-] 

So solubility of solute per 100 g of solute [g g" 1 ] 

t time [s] 

Ug superficial gas velocity [m s" 1 ] 

ua mass fraction of component A in solution [kg A kg" 1 

solution] 

V volume of the liquid in the vessel [m 3 ] 

X impeller clearance from the bottom of tank [m] 

x, y, z exponents in Equation (33) [-] 

Greek letters 

Ap solid-liquid density different [kg m" 3 ] 

a reaction order for cyanide [-] 

reaction order for oxygen [-] 

y reaction order for gold [-] 

6 thickness of the Nernst boundary layer [m] 

Sp thickness of passive layer [m] 

t] reaction order in Equation (29) [-] 

9 exponent for mean diameter of particles in Equation (27) [-] 

H dynamic viscosity [Pa s] 

ju a apparent viscosity according to the Ostwald de Waele 

model [Pa s] 
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density 



[kgm- 3 ] 



Abbreviations 



6FBT flat blade turbine with 6 blades 

B baffled tank 

CBT curved blade turbine 

CBP curved blade paddle 

CCL & Zn solid-liquid separation and zinc precipitation 



CIL 


carbon-in-leach 


CIP 


carbon-in-pulp 


FBP 


flat blade paddle 


ILR 


InLine Leach Reactor 


IS 


inter- stage screens 


MCMC 


Marcov chain Monte Carlo 


MISC 


flotation and gravity concentrate 


ORP 


oxidation-reduction potential 


PBT 


pitched blade turbine 


SCM 


shrinking-core model 


Two-6FBT two stirrers of the type indicated 


UB 


unbaffled tank 
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1 INTRODUCTION 

The chemical element gold, symbol Au, is classified as a noble metal due to its 
inertness to chemical reactions in non-complex media. It does, however, react with 
numerous reagents. It belongs to the same group as copper and silver in the 
periodic table and it is commonly found to be associated with these elements in 
rocks (Juvonen, 1999). Gold is also found in host minerals, typically such as 
calaverite (AuTe2), montbroyite (Au2Te3) and sylvanite (AuAgTe 4 ), in varying 
concentrations and occurs in association with minerals, for example, sulfide and 
copper (Marsden and House, 1992). The average concentration of gold in earth's 
crust is 0.005 g/t, which is much lower than most other metals, for example, 0.07 g 
of silver/t and 50 g of copper/t. The gold content is dependent upon gold minerals 
as well as gold properties, for instance, electrum, specific gravity 16-19.3, is a 
mixture of silver and gold containing 45-75 % gold (Marsden and House, 1992). 
There are many possible methods to recover gold from ores such as leaching, 
gravity concentration and flotation. Leaching by cyanide solutions or gold 
cyanidation, however, has been the main metallurgical process for gold extraction 
for more than one century (de Andrade Lima and Hodouin, 2006). 

Since the invention of the gold cyanidation process in 1887, its chemistry and 
leaching kinetics have been the subjects of considerable investigation and several 
theories have been proposed to explain the reaction mechanism. Various variables 
affecting gold cyanidation, such as dissolved oxygen concentration, free cyanide 
concentration, temperature, pH and particle size, have been studied. Their effects 
on the optimal gold conversion have been investigated since they could result in 
effective improvements at industrial scale, for instance the reduction of the 
operational cost. However, gold cyanidation is a complex system due to the fact 
that gold particles occur as alloys or compounds which are embedded in a mineral 
matrix and galvanic interaction can take place between the phases (de Andrade 
Lima and Hodouin, 2005). A number of these studies remain contradictory in their 
conclusions concerning the mechanism of dissolution. Some studies propose that 
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the diffusion of reactants to the gold surface controls the rate of reaction, while 
others claimed that the chemical reaction is slow (Crundwell et al, 1997). 

For several years mathematical models have been extensively studied in order to 
obtain predictions for the kinetics of leaching reactions. A number of rate 
equations have been published based on different experimental conditions. Some 
of them are mechanistic rate equations and some are empirical ones. The first 
comprehensive study on the variation of the rate of gold dissolution with cyanide 
and oxygen concentration has been described since 1966. However, the limitation 
of the equation in the modeling of industrial gold cyanidation is due to the 
difficulty in measuring the exact surface area of gold in the ore and the Nernst 
boundary layer thickness (Ling et al, 1996). Additionally, numerous researchers 
have claimed that two expressions can be derived from the electrochemical 
mechanism of gold leaching and the shrinking core model with surface passivation 
(Crundwell et al, 1997). The gold leaching process often appears to be operated far 
from an optimum range. For example, excessive reagents (CN~ and O2) 
consumption is expected and the behavior of the process is very dependent on feed 
mineralogy (Ling et al, 1996). Therefore, it is evident that a simple rate equation 
with a minimum number of adjustable parameters is still needed for the modeling 
of gold cyanidation. 

In this work, the aim is to study the effect of leaching parameters on gold 
cyanidation and their limitations. The main purpose is to find the important factors 
which impact the reaction kinetics and to optimize the process of the gold leaching 
by gathering information from existing cyanide leaching processes. The Modest 
(Model Estimation) computer software is used for parameter estimation of the 
mechanistic mathematical models. 
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2 



CHEMISTRY OF GOLD CYANIDATION 



Gold cyanidation has been used as the principle gold extraction technique since the 
late 19 th century. Cyanide is universally used because of its relatively low cost and 
great effectiveness for gold dissolution. Also, despite some concerns over the 
toxicity of cyanide, it can be applied with little risk to health and the environment. 

2.1 Chemistry of cyanide solutions 

Simple cyanide salts, for example, sodium cyanide (NaCN), potassium cyanide 
(KCN) and calcium cyanide [Ca(CN2)], have been widely used as sources of 
cyanide leaching. They dissolve and ionize in water to form their respective metal 
cation and free cyanide ions (CN ~) as presented below: 



Cyanide ions hydrolyze in water to form hydrogen cyanide (HCN) and hydroxyl 
(OH~) ions which also increase pH. At pH about 9.3, half of total cyanide exists as 
hydrogen cyanide and half as free cyanide ions. At higher pH, the total cyanide 
greatly exists as free cyanide ions. Undesirable reactions might occur during 
leaching because hydrogen cyanide, as well as free cyanide, can be oxidized with 
oxygen to form cyanate (CNO~) which does not dissolve gold and thus reduces the 
free cyanide concentration (Marsden and House, 1992). 



NaCN ^ Na + + CN" 



(1) 



CN" + H 2 0 ^ HCN + OH" 



(2) 



4HCN + 3C-2 - 4CNO~ + 2H 2 0 



(3) 



3CN~ + 20 2 + H 2 0 ? 3CNCT + 20H" 



(4) 



2.2 Gold dissolution 



The oxidation of gold is a prerequisite for its dissolution in alkaline cyanide 
solution. Although gold is inert to oxidation, it is widely accepted that, in the 
presence of a suitable complex agent such as cyanide, gold is oxidized and 
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dissolved to form the stable complex ion [Au(CN 2 )~]. Oxygen is reduced and 
hydrogen peroxide is formed as an intermediate product in the first step and 
becomes the oxidizing agent in the second step, leading to the following chemical 
reactions which proceed in parallel (Marsden and House, 1992; Kondos et al, 
1995; Ling et al, 1996; de Andrade Lima and Hodouin, 2005; Senanayake, 2005): 

2Au + 4CN " + 0 2 + 2H 2 0 2Au(CN) 2 ~ + H 2 0 2 + 20H~ (5) 
2Au + 4CN ~ + H 2 0 2 2Au(CN) 2 ~ + 20H~ (6) 

The summation of the two partial reactions is presented in Eq. (7), as proposed by 
Eisner: 

4Au + 8CN " + 0 2 + 2H 2 0 -> 4Au(CN) 2 " + 40H" (7) 

This equation, called Eisner's equation, is stoic hiometrically correct. However, it 
does not completely describe the cathodic reactions associated with the dissolution 
(Marsden and House, 1992). 

The dissolution mechanism has been debated under both acidic and alkaline 
conditions. Dissolution involves an electrochemical process in which the anodic 
reaction is gold oxidation while the cathodic reaction is oxygen reduction. 
Senanayake (2008) has illustrated the gold ion diffusion through the interfaces into 
the solution as shown in Figure 1 . 



metal 


metal-film 
boundary 


film-solution 
boundary 


solution 


Au — 
-e 


film diffusion 

"An** — 


AuCN, 0) 
film growth 


CN" 

► Au(CN) 2 
CN* film dissolution 



Figure 1 Anodic cyanidation model for gold; boundary i: gold-film interface, 
boundary o: film-solution interface (Senanayake, 2008). 
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The important steps during the anodic reaction in the solution phase are presented 
below (Kondos et al, 1995; Ling et al, 1996; Wadsworth et al, 2000): 

(a) Adsorption of cyanide on gold surface 

Au + CN"(aq) i= AuCN " (s) (8) 

(b) Electrochemical extraction of an electron 

AuCN"(s) - AuCN(s) + e" (9) 

(c) Combination of the adsorbed intermediate with another cyanide ion 

AuCN(s) + CN" ^ Au(CN) 2 "(aq) (10) 

where (s) refers to surface adsorbed species and AuCN is a neutral intermediate 
species adsorbed on the surface. 

(d) Overall anodic reaction 

Au(CN) 2 " + e" - Au + 2CN~ (11) 
for which the Nernst equation is: 



E h = -0.60 + 0.1 18 pCN + 0.059 log (C ) V (12) 



The anodic dissolution is accompanied by the cathodic reduction of oxygen 
involving several parallel and series reactions: 

(d) oxygen reduction to H2O2 

0 2 + 2H 2 0 + 2e" - H2O2 E h = 0.69-0.059 pH-0.0295 p0 2 (13) 

(e) hydrogen peroxide decomposition 

2H 2 0 2 - 0 2 + 2H 2 0 (14) 

(f) oxygen reduction to hydroxide ions 

0 2 + 2H 2 0 + 4e" - 40H " E° = 0.401 V (15) 
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The hydrogen peroxide formed is a strong oxidizing agent which may able to 
reduce to hydroxide ions (OH ") as follows: 



However, it has been shown that the reduction of this species is difficult to happen 
and the dissolution rate of gold in oxygen-free solutions containing hydrogen 
peroxide is very slow (Marsden and House, 1992). 

Figure 2 illustrates the major reactions of the two electron processes. The rate 
limiting conditions appear when the diffusion rates of cyanide and oxygen are 
equal. Then the slower diffusion rate of any species will provide the rate limiting 
factor (Marsden and House, 1992). The slow rate has been related to the formation 
of a passive layer on gold surface as well (Zheng et ah, 1995). 



H 2 0 2 +2e" E= 20H 



E° = 



0.88 V 



(16) 



0 2 + 2H 2 0 + 2e - H 2 0 2 + 20H 



2Au + 4CN - 2Au(CN) 2 + 2e 



c 



3 



Cathodic area 



Anodic area 



Solid 



2* 




Au(CN)" 



Figure 2 



Schematic representation of the local corrosion cell at a gold surface 
in contact with an oxygen-containing cyanide solution. ?' a is the 
anodic current; z c is the cathodic current (Marsden and House, 
1992). 
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2.3 Competitive Reactions in Alkaline Cyanide Solution 

Many other elements and minerals are able to dissolve in dilute alkaline cyanide 
solution as well. These competitive reactions may increase reagents consumption 
and consequently reduce the efficiency of gold leaching. The sulfides, for example, 
are dissolved and produce metal cyanide complexes and various sulfur-containing 
species, such as sulfate, sulfide, thiocyanate and thiosulfate ions. For instance, 
pyrrhotite which is the most reactive iron sulfide in alkaline cyanide solution reacts 
with cyanide to form thiocyanate as presented in the following equation 
(Deschenes et al, 1998): 

Fe 7 S 8 + CN" ^ 7FeS + SCN" (17) 

and further oxidation produces Fe(II) cyanide: 

2FeS + 12CN " + 50 2 + 2H 2 0 ^ 2Fe[CN] 6 4 " + 2S0 4 2 " + 40H " (18) 

The Fe[CN]6 4 " complex is stable in the pH range used in cyanidation, as indicated 
in £h-pH diagram for the Fe-S-CN-H 2 0 system at 25 °C (Marsden and House, 
1992). It indicates that the ferrocyanide complex can decompose at higher pH. Iron 
hydroxide is the most likely precipitate after the ferrocyanide complex is 
completely dissociated. The reaction is given as (Rees and van Deventer, 1999): 

2Au + 2 / 3 Fe(CN) 6 4 " + V2O2 + H 2 0 - 2Au(CN) 2 " + 20H " + 2 / 3 Fe 2+ (19) 

Other minerals, such as copper, and metal oxides are soluble to form metal cyanide 
complex in alkaline cyanide solutions as well. Chalcopyrite (CuFeS 2 ) is one of 
common copper minerals showing low solubility in cyanide solution. It is reported 
that concentration of this copper mineral between 300 and 400 ppm in solution 
does not affect gold leaching. All other common copper minerals, for example, 
chalcocite (Cu 2 S) and cuprite (Cu 2 0) have higher solubility. The possibly 
competitive reactions are shown below (Vukcevic, 1997): 

Cu 2 0 + 5CN " + V2O2 + H 2 0 ^ 2Cu(CN ") 2 " + OCN " + 20H " (20) 
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Cu 2 S + 5CN " + V2O2 + H 2 0 s= 2Cu(CN ") 2 " + S + OCN " + 20H " (21) 

These minerals generally consume fewer amounts of cyanide and oxygen than 
sulfide. However, the species formed may affect precipitation reactions and overall 
process efficiency (Marsden and House, 1992). Consequently, several pretreatment 
methods, for example, pre-aeration, pressure oxidation and roasting are applied to 
the leach slurry prior to cyanidation in order to reduce the oxygen and cyanide 
consumption. 



2.4 Oxidation-Reduction Potential (ORP) 



Redox potential or ORP can be used to explain the stabilities of metals and other 
species in aqueous solutions. ORP is related to the potential-pH diagrams, also 
called £h-pH or Pourbaix diagrams, by Nernst equation. Each line on the £h-pH 
diagram represents the condition where the activities of reactants and products of 
the considered reaction are in equilibrium. Figure 3 illustrates the £h-pH conditions 
applied in industrial gold extraction process. This Pourbaix diagram indicates that 
the reaction of the Au(I) complex takes place more readily than the Au(III) 
complex since the Eh value of the Au(I) complex reaction is more negative than the 
other. Wadsworth et al. (2000) presented that gold dissolution in alkaline cyanide 
solution (Au(CN)2~ ) was found to occur at the redox potential range of -0.4 and 
-0.7 V. Figure 3 also indicates that the potential difference between the lines of 
gold oxidation and oxygen reduction reactions is maximized at pH value of about 
9.5. The Pourbaix diagram is also able to explain the role of other minerals 
competing with gold cyanidation. 
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Figure 3 Potential-pH diagram for the system AU-H2O-CN" at 25 °C. 

Concentrations of all soluble gold species = 10" 4 M (Marsden and 
House, 1992). 

Figure 4 illustrates the iron mineral system, one of the important competitve 
minerals in alkaline cyanide solution. Figure 4 shows that the iron(II) complex can 
be formed readily. According to the comparison between the Figures 3 and 4, it 
can be seen that the Fe(II) complex is more reactive in cyanide solution than Au(I) 
at the same value of pH since the E h value of the Fe(II) complex is more negative. 
Therfore, the iron minerals might be reacted with cyanide thus retarding gold 
dissolution. However, the E h value also depends on process conditions. 
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Figure 4 £h-pH diagram for the Fe-S-CN-H 2 0 system at 25°C (Marsden and 
House, 1992). 
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3 PHENOMENA IN GOLD LEACHING 

The most important reactions in hydrometallurgical gold extraction processes are 
heterogeneous, involving the transfer of metals and minerals between solid and 
liquid phases. Heterogeneous reactions are controlled either by the inherent 
chemical reaction kinetics or by the rate of mass transfer of the individual reacting 
species across a phase boundary (Marsden and House, 1992). 

3.1 Reaction Kinetics 

Many researchers have attempted to model the kinetics of gold dissolution. The 
rotating gold discs have been used in fundamental studies to determine the rates at 
a constant surface area, assuming that the surface roughness does not change 
during the course of reaction. It was concluded that the rate of pure gold 
dissolution relies on the rate of film diffusion of cyanide ions or dissolved oxygen 
towards the gold surface as shown below (Kondos et al, 1995; Ling et al, 1996): 

-m±~ 2D ™- D °^- C °. molm - v (22) 

Ad, Si(D c „-C cll _) + (4D 0i C O: )} 

where A surface area of gold disc in contact with aqueous phase, m 2 
C CA ,_ cyanide concentration, mol/m L 3 



oxygen concentration, mol/mi/ 



D CN _ diffusion coefficients of cyanide, m 2 /s 



2 

D 0 diffusion coefficient of oxygen, m /s 



Nau amount of gold in particles, mol 
t time, s 

6 The Nernst boundary layer thickness, m 
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From Equation (22), it can be seen easily that when D CN - C CA ,_ <4.D 0o C 0i or low 

cyanide concentration, gold dissolution rate depends primarily on that cyanide 
concentration: 

Au 1 ^CN „ , „ , -2 -1 / 0 o N 

-^T = r^ c »- =i » c »- molm s (23) 

where k CN _ rate constant, s" 1 

Similarly, for high cyanide concentration, the gold dissolution rate becomes mainly 
dependent on the oxygen concentration: 

- = 2^- C 0 = k 0 C 0 mol m" 2 s 1 (24) 

Adt 8 2 2 2 

where /: 0 rate constant, s" 1 

In practice, a high level of cyanide has been maintained rather than a high 
dissolved oxygen level in the solution. Consequently the majority of the mills 
operate at cyanide levels such that gold dissolution is dependent on the dissolved 
oxygen level (Kondos et al, 1995). However, the use of Equations (22)-(24) is 
limited in modeling an industrial process due to the difficulty in measuring the 
exact gold surface area in the ore and Nernst boundary layer thickness (Ling et al, 
1996). 



An expression of second order rate equation was proposed by Nicol et al. (1984), 
as presented in Equation (25). This equation, the so called Mintek equation, was 
presented as an empirical equation based on the leaching behavior of several South 
African gold ores and has the following form (Ling et al, 1996; Crundwell and 
Gordorr, 1997): 

dC Au _ , (n ^~ N 2 ,wl ,-1 



dt 



= k(C Au -CZY mgkgV (25) 



where k overall rate constant, kg/mg s 
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C Au gold concentration in particles, mg/kg 

C°° Au gold concentration in particles after infinite leaching time, mg/kg 

This equation neither accounts for the reagent concentration, C CA ,_ and C Qi , nor for 

the particle size (Ling et ah, 1996; de Andrade Lima and Hodouin, 2005). 
Consequently, many researchers have further endeavored to present the kinetics 
involved a term of the reactants. 

Ling et al. (1996) proposed a single rate equation on the assumption of pseudo- 
homogeneous ore behavior. This rate equation can be expressed in terms of mass 
per mass of ore instead of mass per surface area of gold units, as follows: 

-^ L = kC c N a C o/(C Au -CZY mgkg 1 h- 1 (26) 

where C Au gold concentration in particles, mg/kg 

C Au gold concentration in particles after infinite leaching time, mg/kg 
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C CA ,_ cyanide concentration, mg/dm 

C 0i oxygen concentration, mg/dm 3 

k overall rate constant, unit depends on reaction orders 
a reaction order for cyanide 
reaction order for oxygen 
y reaction order for gold 

In the model, a reaction order of 1.5 for gold was found to be the best fit to 
experimental data. They also found that the overall rate constant, k, equals to 
0.0016 + 0.0002, a equals to 0.81 + 0.10 and /J equals to 0.73 + 0.09 within the 
examined range of C CA ,_ 25-156 mg/dm 3 , C 0i 8.5-40 mg/dm 3 , Ca u 4.7-13.2 mg/kg 

and pH 10. This equation does not account for the ore particle size. However, it 
was claimed to be used for preliminary process analysis. 
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de Andrade Lima and Hodouin (2005) has recently developed the rates equation 
describing the gold dissolution, using the general pseudo-homogeneous kinetic 
equation, as presented in equation (26), considering also the particle size. The 
overall rate constant, k, is a function of the average diameter of the ore particles, as 
follows: 

k = c l -c 2 d e (27) 

where c l , c 2 coefficients, units depend on reaction order 

d average diameter of particles, um 

9 exponent for average diameter of ore particles 

The gold dissolution kinetic equation shows that the reaction orders are 
approximately two for gold, one for free cyanide and a quarter for dissolved 
oxygen. The constant is a decreasing cubic function of the average particle 
diameter and may be controlled by the particle volume. The rate equation for gold 
dissolution is shown below: 

d c 

_^ = (L13x1 0-3_4.37 X io- 11 j 2 - 93 )x (C A „, e -C2, c ) 213 

dt 

^ 0.961 „ 0.228 , -1 , -1 /00 , 

xC CN - C 0i mgkg h (28) 

This equation was obtained by using free cyanide concentrations 260 and 650 
mg/dm 3 , dissolved oxygen concentrations 8 and 40 mg/dm 3 , gold content in the ore 
1.5-2.3 mg/kg and pH 12. This equation was proved to fit well the fine size 
fractions because gold particle segregation was negligible and small samples could 
adequately represent the system, while this is acceptable for the coarse particles. 
Since cyanide is not a selective leaching agent, it can react with other chemical 
species, such as zinc, iron, silver and copper minerals. Thus, de Andrade Lima and 
Hodouin (2005) also investigated the kinetic model for the cyanide consumption as 
a function of cyanide concentration and particle size. The following pseudo- 
homogeneous model describes the consumption process: 
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- r ar= k ar C ar n mg dm 3 h 1 (29) 



(30) 



CN ~ d 9 -c 4 
where C CA ,_ cyanide concentration, mg/dm 3 

c 3 , c 4 coefficients, units depend on reaction order 
d average diameter of particles, um 

k overall rate constant which is a function of the average diameter of 

particle, (dnrVmg) 271 h" 1 
9 exponent for average diameter of ore particles 

The estimated cyanide consumption kinetic shows that the reaction order is 
approximately three for free cyanide. The rate constant is a decreasing reciprocal 
square-root function of the particle diameter. The free cyanidation model is then 
given as follows: 



r CN- 



( \ 

1.69xl0" 8 
v J 0547 -6.40 y 



<V 371 mgdm^h- 1 (31) 



It was also found that the cyanide consumption kinetics are faster for the small 
particles due to the liberation of the cyanide-consuming species found in the ore 
(de Andrade Lima and Hodouin, 2005). 



The electrochemical reaction mechanism was also used to describe the rate of gold 
leaching by Crundwell and Godorr (1997) as given below: 



dN - =Uk a rr(k c C o r molm-V 1 (32) 



Adt F 



where F the Faraday constant, 96,485.34 C/mol 

k a rate constant of anodic half-reaction, C m/mol s 
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rate constant of cathodic half-reaction, C m/mol s 



The rate expression which is derived from the anodic and cathodic half-reactions is 
one-half order in the concentrations of cyanide and oxygen. The cyanide 
concentration determines the rate of anodic dissolution of gold while the oxygen 
concentration determines the cathodic reduction of oxygen (Crundwell and Godorr, 
1997). 

3.2 Gas-Liquid Mass Transfer 

Mass transfer from gas to liquid phase has a significant role in cyanidation using 
oxygen in the stirred tank reactor. Gas is introduced in the form of bubbles into the 
liquid by an appropriate distributor. The gas-liquid mass transfer is normally 
described by the volumetric mass transfer coefficient, kifl s" 1 where k L is the mass 

2 3 

transfer coefficient, m/s, and a is the specific interfacial area, m /m . The most 
important characteristics affecting the gas-liquid mass transfer are the energy 
dissipation, the gas hold-up and the bubble size. These variables are a function of 
the vessel geometry, mainly the stirrer and gas distributor, and operational 
conditions such as power input and gas flow. Physical properties of solution and 
gas phase, for example, viscosity, surface tension and density also have effect on 
the parameters of gas-liquid mass transfer (Garcia-Ochoa and Gomez, 2004). 

Researchers have attempted to present the volumetric mass transfer coefficient by 
using the experimental data and empirical correlations obtained in stirred tank 
reactor. Garcia-Ochoa and Gomez (2004) have presented a correlation for the kifi 
as follows: 



k L a 



= mU x G (PIV) y n z 



s 



-l 



(33) 



where m 



constant 



P 



power input under gassed conditions, W 



U G 



superficial gas velocity, m/s 
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V volume of liquid in vessel, 
ju viscosity, Pa s 

Other equations proposed substitute the average power input per volume, P/V, by 
the effect of stirrer speed, N. Garcia-Ochoa and Gomez (2004) also have presented 
the volumetric mass transfer coefficient for different stirred tank volumes and for 
both Newtonian liquid, which is air-water system, and non-Newtonian liquid, such 
as air-sodium sulfite solution system. The exponents (x, y and z) are shown in 
many different dimensional correlations based on tank configurations, properties of 
liquid and operational conditions as given in Table 1. 

The stirred tanks presented are single and multi impeller systems such as flat blade 
turbine, flat blade paddle, pitched blade paddle and two stirrers of flat blade turbine 
with four blades. Garcia-Ochoa and Gomez (2004) have found an increase in kha 
with the increasing of superficial gas velocity, power input (stirred speed) and 
liquid viscosity. Fan and Herz (2007) have proposed that the kia values between 
0.05-0.5 s" 1 can be achieved in most well-designed stirred tank reactors. They have 
found that the k L a value for C02-Na2C03 solution system was 0.06 s" , where the fcz, 
and a values were 7.2 xlO" 4 m/s and 84 m 2 /m 3 , respectively. Lorenzen and 
Kleingeld (2000) have presented the kia value for carbon dioxide gas-sodium 
hydroxide (C02-NaOH) solution system in the high intensity gas-liquid jet reactor. 
The results have been shown to exhibit high mass transfer coefficient, interfacial 
areas and volumetric mass transfer coefficient compared to the conventional 
system. The k L value was up to 3xl0~ 3 m/s while the value of a was ranging from 
2000 to 16000 m 2 /m 3 . Consequently, the kia was up to 25 s" 1 which is higher than 
the kha for conventional processes in the region of 0.025-1.25 s" 1 . Lorenzen and 
Kleingeld (2000) have claimed that these results can be compared to the possibly 
enhanced mass transfer rates when using oxygen or air during cyanidation. 
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3.3 Solid-Liquid Mass Transfer 

Solid-liquid dispersion in mechanically agitated vessels or stirred tank reactors has 
been widely studied. The main reason for applying the mechanical agitation is to 
ensure that all of the surface area available for mass transfer is utilized. Evaluation 
of the solid-liquid mass transfer coefficient, k S u has often been performed by solid 
dissolution. This coefficient is dependent upon the homogeneity, which is the 
function of geometric configurations (type of impeller, impeller size to tank size 
ratio or dWdi, impeller location or X/dj), operating parameters (impeller speed, 
power input, particle loading) and physical properties of the particles and fluid 
(viscosity, solid-liquid density difference, particle size and shape). The turbulence 
intensity close to the impeller is the highest and decreases as the distance from the 
impeller increases (Pangarkar et al, 2002). 

Pangarkar et al. (2002) have presented a correlation for the solid-liquid mass 
transfer coefficient in three-phase (solid-liquid-gas) system, where the type and the 
location of the sparger and gas flow rate are taken into account. The correlations 
are presented on the basis of classified approaches which are dimensional analysis, 
Kolmogoroff s theory of isotropic turbulence, slip velocity theory and analogy 
between momentum transfer and mass transfer. The dimensional approach-based 
correlations are presented in terms of Sherwood number based on tank diameter, 
Shi (ksLdj/D m ), and particle diameter, Shp(kshdp/D m ). The correlations are given in 
Table 2. Pangarkar et al. (2002) suggested a correlation with many variables such 
as liquid viscosity and density, impeller speed, dr/dj and particle size. Some 
correlations are presented only on the terms of Reynolds number based on impeller 
diameter, Re^ (Ndvplp.) an d tank diameter, Rej (Ndj 2 p/p) and Schmidt number, Sc 
Qj/pD m ), while some have several variables. The constants and the exponent of the 
Reynolds number and Schmidt number vary with geometric configurations and 
operating parameters. Thus, the application of the correlation based on dimensional 
approach is limited since it predicts widely different &sl values and can only be 
applied to geometrically similar configuration under similar range of operating 
parameters. The estimation of k S L values from momentum transfer-mass transfer 
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approach-based correlations, which correlate relative particle suspension, is 
claimed to be reliable and simple. In this case, k S L is a function of impeller speed 
and Schmidt number. Other approaches, to correlate ksu are also presented as a 
function of Reynolds number and Schmidt number with different exponent values. 
However, there are also limitations in the Kolmogoroff's theory about geometric 
variables and in the density difference and slip velocity theory about turbulence 
characteristics (Pangarkar et al, 2002). 

In the case of fine particles (dp < 100 um), which are commonly found in processes 
such as crushed ores in leaching tank, the mass transfer coefficient increases 
greatly with decreasing particle size. The ksh values are claimed to be independent 
of the density difference and tank diameter and can be correlated as Sh P =f(Re £ , Sc) 
based on the Kolmogoroff's theory. In the correlation, the Sherwood number based 
on the particle diameter, Sh P , is a function the Reynolds number and Schmidt 
number (Pangarkar et al, 2002). 
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Table 2 Correlations based on dimensional approach proposed by different 

authors (Pangarkar et al, 2002). 
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4 VARIABLES AFFECTING THE PHENOMENA 

Leaching parameters have been widely studied to optimize the performance of the 
cyanidation process. Performance of process depends on variables such as the 
concentrations of dissolved oxygen, free cyanide and compounds that react with 
cyanide, pH, particle size and operating conditions. These parameters affect on the 
gold dissolution rate, gold extraction and cyanide consumption. Many studies have 
been carried out under different experimental conditions in order to investigate the 
effect of leaching parameters and obtain the optimum process conditions. Some of 
those are listed below. 

Deschenes and Wallingford (1995) experimented with a sulfide bearing gold ore 
containing 6.75 g/t Au. The ore was grinded to the size of 74 urn. The stirring 
speed was kept constant at 400 min" 1 . The ore was pre-leached with air for four 
hours followed by cyanidation. 95.8% Gold recovery was obtained under condition 
of free cyanide concentration 480 ppm, lead nitrate 200 g/t, pH 10.5 and 16 ppm 
O2. The cyanide consumption in 24 hours was 0.51 kg/t 

Ling et al. (1996) presented that the maximum gold dissolution achieved was about 
93 % under ambient temperature 20 °C and air sparging at pH 11, dissolved 
oxygen of 8.5 ppm and cyanide concentration of 104 ppm. The agitation speed was 
maintained at 750 min" 1 for 24 hours. Ore sample in this case was from a northern 
Ontario mining site and had a mean size of 37 um. 

Deschenes et al. (2003) performed experiments with the dissolved oxygen 
concentration at 10 ppm and pH at 11.2. A sample of high grade ore 77.8 g/t 
received from Goldcorp Red Lake Mine, Canada was milled to the size of 37 um. 
The cyanide consumption for a 24-h leach was 0.55 kg/t using 300 ppm NaCN. 
The experiment was carried out at a constant rate of 400 min" 1 . It was found that 
adding 100 grams of lead nitrate per ton to the pre-leach circuit increases gold 
recovery. At these conditions, about 96% gold extraction was achieved. 
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Ellis and Senanayake (2004) investigated gold leaching in multi-mix system at 
ambient temperature, pH about 10, NaCN about 200 mg/dm 3 and dissolved oxygen 
about 20 mg/dm 3 . In this study, a pyrrhotite-rich ore, 63 um, from the Bounty Gold 
mine located in Western Australia was used. They concluded that pre-oxidation for 
12 h followed by cyanidation of C CA ,_ I C Qj molar ratio at about 12 provided the 

best gold extraction in the first cyanidation tank. The percentage of Au extracted in 
the first tank reached a maximum of 83% when the product C x C Q was about 

5.6 mmol 2 /(dm 3 ) 2 . 

Consolidated Murchison located in South Africa is one of the gold industries using 
cyanide leaching together with oxygen. In the process, the ore is crushed to the size 
of about 75 um. Around 1983, 77% gold recovery was achieved under conditions 
at pH 7 by lime, 2% NaCN, 52% solids and ambient temperature of 22-38°C. 
Leaching is performed in a pipe reactor at a pressure of 9 MPa, with pure oxygen 
injected into the reactor at 12 MPa. The reactor provided a retention time of 15 
min. Cyanide and oxygen consumptions were 57 and 20 kg/t, respectively 
(Marsden and House, 1992). 

Table 3 summarizes the optimum operating conditions from various experiments 
as written above. 
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4.1 Concentration of Oxygen 

Oxygen is one of the key reagents in conventional gold cyanidation as presented in 
Equation (7). The concentration of oxygen determines the cathodic reduction of 
oxygen. It may be supplied to the system by air, pure oxygen or as enriched air. 
The oxygen concentration for processes using air is determined by the conditions 
of temperature and pressure that the process operates under. The saturated 
dissolved oxygen concentration is about 8.2 mg/dm at 25°C (Marsden and House, 
1992). Ellis and Senanayake (2004) have suggested that oxygen injection, instead 
of air, provides a high oxygen concentration in gas to increase the oxygen 
solubility and can achieve high overall gold extraction. However, Kondos et al. 
(1995) found that the oxygen consumption is markedly higher with the use of 
oxygen gas than air. A large number of small oxygen bubbles dispersed in the 
slurry long enough and deep enough gives adequate oxygen concentration for gold 
dissolution (Ellis and Senanayake, 2004). Knowing the oxygen demand the process 
can be optimized in order to obtain higher gold recovery, shorter retention time, 
low cyanide consumption which reduces environmental concerns relating to 
cyanide in the effluents. However, different mineralogical compositions need 
different amount of oxygen, for example, pyrite-rich ore has a higher oxygen 
demand than ore containing albite, quartz and chlorite as the main minerals 
(Kondos etal, 1995). 

Many researchers have concluded that increasing the dissolved oxygen 
concentration increases the rate of dissolution (Marsden and house, 1992; 
Deschenes et al, 2003; Ellis and Senanayake, 2004). Deschenes et al. (2003) have 
presented the effect of dissolved oxygen on gold extraction as illustrated in Figure 
5. They discussed that dissolved oxygen concentration has no significant effect on 
cyanide consumption and further discussed that faster leaching kinetics were 
observed by using higher dissolved oxygen concentrations. 
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Time (h) 

Figure 5 Effect of oxygen addition on cyanidation: pH 11.2, 500 ppm NaCN, 
24 h (Deschenes et al, 2003). 



However, Ling et al. (1996) have presented Figure 6 and concluded that high 
dissolved oxygen level did not have significant effect on gold dissolution rate at 
the cyanide concentration level of 104 ppm. 




Time<h) 

Figure 6 Effect of dissolved oxygen concentration on gold conversion (Ling 
et al, 1996). 

4.1.1 Solubility of Oxygen 

The oxygen solubility depends on many factors such as temperature, pressure and 
ionic strength (Ellis and Senanayake, 2004). Ling et al. (1996) have discussed that 
the oxygen level in the solution is mainly affected by temperature and pressure. 
The solubility of gases in aqueous solutions generally increases with decreasing 
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temperature and increasing gas partial pressure. The ionic strength is relative low 
and its effect on gas solubility is not important (Ling et al, 1996). 

4.2 Cyanide Concentration 

Many studies have presented that gold extraction increases with increasing cyanide 
concentration (Marsden and House, 1992; Kondos et al, 1995; Ling et al, 1996; 
Wadsworth et al., 2000; Deschenes et al., 2003). Figure 7 shows the effect of 
cyanide concentration on gold conversion. 



Ling et al. (1996) have found that the gold dissolution rate was not very sensitive 
to the change of cyanide concentration at high cyanide levels. This was also found 
by Deschenes et al. (2003) who presented in Figure 8 that the use of 400-500 ppm 
NaCN resulted in similar overall gold extraction. 
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Figure 7 



Effect of cyanide concentration on gold conversion (Ling et al., 
1996). 
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Figure 8 



Effect of cyanide concentration on leaching. Pre-leaching: pH 11.2, 
8 ppm 0 2 , 100 g/t Pb(N0 3 ) 2 , 12 h; cyanide: pH 11.2, 10 ppm 0 2 
(Deschenes et al, 2003). 



Ellis and Senanayake (2004) have proposed that the rate of gold leaching increases 
with increasing concentration of cyanide but becomes independent of cyanide 
concentration when it exceeds 0.075% KCN (or 0.06% NaCN, equivalent to 600 
ppm). Kondos et al. (1995) have concluded that using an excess of cyanide results 
in unnecessary cyanide consumption and has no beneficial effect on gold 
extraction. The excess cyanide will consume more cyanide because of the 
formation of cyanocomplexes from impurities. However, a high concentration of 
cyanide may be required due to the competition of other species such as ore 
containing soluble sulfides (Mardens and House, 1992). Furthermore, the decrease 
in cyanide concentration decreases cyanide consumption as well. The use of lower 
cyanide level can reduce the effluent treatment cost at the plant scale (Ling et al, 
1996; Deschenes et al, 2003). Therefore, lower or higher cyanide concentrations 
than the optimum value can provide negative consequence to the gold leaching 
process. 



The addition of lime (CaO) is typically used to maintain pH in gold leaching 
circuit (Ling et al., 1996; Ellis and Senanayake, 2004). The gold dissolution rate is 
expected to reduce with increasing pH since the adsorption of OH ~ ion onto gold 
surface decreases the surface available for cyanide leaching (Kondos et al, 1995; 



4.3 pH 
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Ling et al., 1996). This is also due to an increase in the rate of competing reactions 
such as dissolution of sulfides and other reactive species (Marsden and House, 
1992). Parga et al. (2007) have presented the effect of pH on gold extraction as 
given in Figure 9. It can be seen that the gold extraction for high pressure system is 
reduced when pH increases as well (Parga et al, 2007). 




Figure 9 The effect of pH on gold extraction. Condition: 20% solid, 0.6 MPa, 
300 min" 1 , 80°C, 1% NaCN, 1 h (Parga et al, 2007). 



Ellis and Senanayake (2004) have reported that cyanide leaching at high pH is able 
to minimize the formation of hydrogen cyanide causing cyanide loss as HCN gas 
as presented in Equations (2)-(4). A high pH level also results in lower cyanide 
consumption as presented in Figure 10 (Ling et al, 1996). However, the optimum 
pH for leaching should be derived independently for each ore type and leaching 
system (Marsden and House, 1992). 
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Figure 10 Effect of pH on cyanide consumption (Ling et al, 1996). 
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4.4 Particle Size 



Ling et al. (1996) have studied the effect of particle size on gold conversion and 
found that smaller particle size can improve gold dissolution rate. This is because 
smaller particle size provides larger contacting surface area between solid and 
liquid and thus increases the rate. Similar results have been obtained by de 
Andrade Lima and Hodouin (2005) who presented the residual gold concentration 
in ore particle as a function of the ore particle size, as shown in Figure 1 1. 



04 




Average Size (urn) 

Figure 1 1 The residual gold concentration as a function of the ore particle size 
(de Andrade Lima and Hodouin, 2005). 



de Andrade Lima and Hodouin (2005) have further presented the gold dissolution 
rate for different ore particle sizes, as displayed in Figure 12. However, Kondos et 
al. (1995) reported that increasing the fraction of finer grinding will also increase 
cyanide consumption and does not improve gold extraction, probably because 
residual gold is present in fine occlusions or in solid solution in the mineral matrix. 
Marsden and House (1992) also presented that the rate of dissolution may decrease 
with decreasing particle size due to the increased rate of competing and reagent- 
consuming side reactions. Moreover, the grinding costs also increase with the 
decrease in particle size (Ellis and Senanayake, 2004). 
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Figure 12 Gold dissolution rate for ore average particle size 30 and 177 um: 
(a) as a function of the dissolved oxygen concentration; (b) as a 
function of the free cyanide concentration (de Andrade Lima, 2005). 



4.5 Temperature 



Temperature has a significant effect on gold dissolution rate and oxygen solubility. 
Higher temperatures result in higher reaction rates of the metals in the solvent. The 
viscosity of liquid is lower and the diffusivities larger at higher temperature. The 
gold dissolution is, therefore, aimed to increase with increasing temperature as 
present in Figure 13 by Marsden and House (1992). Figure 13 indicates that the 
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dissolution rate of gold in cyanide solution reaches maximum at approximately 
85°C. In contrast, the gas solubility tends to decrease with increasing temperature. 
Furthermore, temperature which is too high leads to leaching of excessive 
undesirable minerals as well. Consequently, the ambient temperature is usually 
applied (Marsden and House, 1992). However, elevated temperatures can be 
applied with high pressure in order to improve gold extraction as discussed in 
Chapter 5.3. 
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Figure 13 Effect of temperature on gold dissolution rate in aerated 0.25% 
KCN solution (Marsden and House, 1992). 

4.6 Pressure 

Pressure can improve gold dissolution in terms of gas solubility. Increasing 
pressure increases the solubility of gas (Ling et al., 1996; Ellis and Senanayake, 
2004). High pressure, together with high temperature, is claimed to improve gold 
extraction as discussed in Chapter 5.3. 

4.7 Slurry Density 

Leaching is normally performed at slurry densities between 35% and 50% of 
solids. This is dependent upon such factors as the specific gravity of solids, particle 
size and the presence of minerals affecting slurry viscosity, such as clays. Mass 
transfer phenomena are maximized when slurry density is low. However, an 
increase in slurry density can increase the solids retention time in the reactor. 
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Furthermore, the high slurry density can reduce the reagent consumption since the 
slurry concentration can be obtained by smaller volume of solution per unit mass 
of material (Marsden and House, 1992). 

4.8 Mixing 

Agitation is one of the important factors in gold leaching since gold dissolution 
rate depends on the diffusion layer thickness and mixing characteristics of the bulk 
solution (Marsden and House, 1992). Therefore, agitation should be sufficient to 
suspend all the particles in the slurry. Increasing the agitation speed can enhance 
the dissolution rate since the intense mixing reduces diffusion layer thickness thus 
substantially improves the mass transfer rates of cyanide and oxygen and allows 
possible saturation of slurry with oxygen to achieve high dissolved oxygen levels 
(Ellis and Senanayake, 2004). Figure 14 presents the effect of impeller speed on 
the dissolution rate of gold disc at different oxygen concentrations. 
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Figure 14 Effect of agitation on the dissolution rate of gold disc (Marsden and 
House, 1992). 
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4.9 Presence of Sulfide minerals and Other Ions in Solution 

Many metals and minerals have a retarding effect on gold leaching rates. These 
metals and minerals are also able to dissolve in alkaline cyanide solution and 
consequently consume cyanide and oxygen. However, it was found that some ions 
can enhance the gold dissolution rate (Marsden and House, 1992; Aghamirian and 
Yen, 2005; Liu and Yen, 1995). 

4.9.1 Effect of Sulfide Minerals 

A sulfide mineral contains sulfide (S 2 ) as the main anion. Most metal sulfides 
decompose readily in aerated alkaline cyanide solution. Sulfide minerals are also 
oxidized to form thiocyanate or to react with oxygen, forming sulfide and sulfate 
(Marsden and House, 1992). Sulfide ions are strongly adsorbed onto gold surfaces 
forming gold sulfide film which passivates the gold surface, and hence retards gold 
dissolution (Kondos et al, 1995) However, concerning the effect of sulfide 
minerals on gold dissolution, it has been found that it is dependent upon mineral 
types and oxygen levels in the solution (Liu and Yen, 1994). 

4.9.1.1 Effect of Pyrite, Chalcopyrite and Pyrrhotite 

Pyrite (FeS2) and Pyrrhotite (Fe(i_ x )S) can introduce various sulfur species and 
ferrous cyanide ions into the solution, while chalcopyrite forms various cuprous 
cyanide compounds (Aghamirian and Yen, 2005). It has been concluded by many 
researchers that sulfide ions have negative effect on gold dissolution causing the 
surface of gold to be passivated (Marsden and House, 1992; Liu and Yen, 1994; 
Aghamirian and Yen, 2005; Dai and Jeffrey, 2006). The effects of sulfide minerals 
on gold dissolution are present in Figure 15. 
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Figure 15 Effect of pyrite, chalcopyrite, pentlandite and oxygen on gold 
dissolution. Condition: 0.25 g of NaCN/1, 400 min" 1 , pH 10.5, 22°C 
(Liu and Yen, 1994). 

Liu and Yen (1994) have concluded that gold dissolution rate is accelerated when 
oxygen concentration is increased in the presence of chalcopyrite and pyrite. 

4.9.1.2 Effect of Galena and Arsenopyrite 

The effect of galena (PbS) on gold dissolution rate is quite different from those 
sulfide minerals discussed earlier. The addition of galena causes the gold 
dissolution to be increased under atmospheric condition at 8 ppm O2. This is 
because a small amount of metal ions is dissolved. The positive effect is also found 
in Arsenopyrite (FeAsS). Contrastingly, the gold dissolution rate reduces under 
oxygen-enriched condition, 20 ppm O2, where large amounts of galena and 
arsenopyrite are dissolved. The rate of gold dissolution is, therefore, reduced. The 
rate is found to be almost zero in galena slurry at 32 ppm O2 (Liu and Yen, 1994). 
Figure 16 illustrates the effect of galena and arsenopyrite on gold dissolution. 
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Figure 16 Effect of arsenopyrite, galena and oxygen on gold dissolution (Liu 
and Yen, 1994). 



It might be concluded that the accelerating and retarding effects of sulfide minerals 
are associated with the solubility of minerals in cyanide solution. The cyanide 
consumption is increased with increasing solubility of minerals and thus the gold 
dissolution rate decreases. 



4.9.2 Effect of Copper Ions 

Copper minerals can dissolve to form a variety of copper cyanide complexes such 
as CuCN, Cu(CN) 2 \ Cu(CN) 3 2 " and Cu(CN) 4 3 " at appropriate pH values (de 
Andrade Lima and Hodouin, 2005). The distribution of copper ions and 
compounds depends on the ratio of cyanide to copper and pH. The formation of 
Cu(CN) 2 " is preferred at low cyanide concentration and at pH below 6, while high 
cyanide concentration and high pH are related to Cu(CN) 3 2 ". Hence, copper 
dissolution is undesirable since it consumes cyanide and oxygen (Marsden and 
House, 1992). An adequate free cyanide concentration in the solution can prevent 
the effect of copper. The ratio of cyanide to copper, approximately 3, is presented 
to be sufficient to remain high gold leaching rate (Marsden and House, 1992; Dai 
and Jeffrey, 2006). When the cyanide to copper ratio at 2.5 is found, the gold 
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leaching rate is low (Dai and Jeffery, 2006). Aghamirian and Yen (2005) also 
presented that copper cyanide at moderate concentrations has a slightly negative effect 
on gold dissolution. 

4.9.3 Effect of Iron Ions 

Iron oxides, such as hematite (Fe 2 03) and magnetite (Fe 3 04) and iron silicates, are 
normally insoluble in alkaline cyanide solutions. Iron carbonates, for instance 
siderite (FeCOs), and other complex carbonate minerals can decompose in low 
alkalinity solution or pH below 10 but are unreactive at high pH values. Iron 
sulfides greatly decompose in alkaline cyanide solutions to form iron cyanide 
complexes and various sulfur species. They are also passivated by the formation of 
a layer of iron (III) hydroxide on the surface. Consequently, gold particles might be 
coated and leaching efficiency may be reduced (Marsden and House, 1992). Iron 
sulfide, especially pyrrhotite (Fei- x S), is a large consumer of cyanide and oxygen 
since it may form the very stable ferrocyanide ion [Fe(CN)6] 4 " (de Andrade Lima 
and Hodouin, 2005). However, it has been indicated that ferrocyanide ions at high 
concentration have a positive effect on gold dissolution (Aghamirian and Yen, 
2005). 

4.9.4 Effect of Lead Ions 

In Several studies it has been presented that addition of lead solution for example 
lead nitrate, Pb(NC>3)2, is able to enhance gold leaching kinetics. Deschenes and 
Prud'homme (1997) have presented that lead nitrate can act in different ways. 
They observed that lead nitrate can activate the surface of a passivated gold 
particle, prevent a passivation film formation on the gold surface and also 
precipitate the soluble sulfides. Kondos et al. (1995) have presented that sulfide 
mineral, which is distributed in the solution, will be precipitated by excess lead as 
PbS. Lead ion may also react with gold sulfide film making the passive gold 
particle to become active again in the absence of sulfide formation. Deschenes et 
al. (2000) have also demonstrated the effect of lead ions in the form of lead nitrate, 
lead sulfide and lead sulfite, which are reacted with gold and accelerate gold 
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dissolution. Figure 17 illustrates the effect of adding lead nitrate on gold ore with 
pyrite, pyrrhotite and chalcopyrite. 




Time (h) 

Figure 17 Effect of lead nitrate addition on the leaching of a gold ore with 
pyrite, pyrrhotite and chalcopyrite (Deschenes et al. 2000). 

Recently, Breuer et al. (2008) have also found that sulfide ion concentration is 
decreased fast with the addition of lead ions, as shown in Figure 18. They further 
concluded that lead sulfide might act as a catalyst for sulfide oxidation releasing 
the lead ions to further react with sulfide ions. Deschenes et al. (2003) have shown 
that the addition of lead nitrate does not affect cyanide consumption. 
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Figure 18 Effect of lead (II) on sulfide oxidation in the absence and presence 
of cyanide (20.4 mM NaCN, 10 mM S ", oxygen sparged) (Breuer 
et al, 2008). 
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However, Kondos et al. (1995) also presented that higher concentrations of lead 
nitrate can reduce gold extraction as shown in the Figure 19. They further 
discussed that the addition of lead nitrate together with oxygen provides the most 
effective gold extraction. 
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Figure 19 Effect of oxygen and lead nitrate on gold leaching. Condition: 
[CN - ] = 730 ppm, pH =11.5, 80% - 74 urn (Kondos et al, 1995). 

4.10 Presence of Activated Carbon in Pulp 

The activated carbon is used as adsorbent in CIP and CIL processes as described in 
Chapter 5. It may enhance gold leaching from ores by continually removing the 
leached species from the solution, redistributing the solid/liquid equilibrium and 
favoring further gold dissolution (Rees and van Deventer, 2000). The chemical and 
physical properties of the activated carbon affect the adsorption kinetics and the 
equilibrium gold-loading capacity. According to Ibragimova et al. (2007), the most 
important characteristics of carbon to be predetermined are the kinetic parameter of 
gold sorption on carbons and the mechanical strength. The adsorption ability has 
been shown to be dependent upon temperature and cyanide concentration. Marsden 
and House (1992) have presented the effect of these parameters, as shown in Table 
4. It has been observed that increasing the temperature and cyanide concentration 
reduces the capacity of gold adsorption on carbon. 
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Table 4 Effect of temperature and sodium cyanide concentration on gold 

loading (Marsden and House, 1992). 
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Rees and van Deventer (2000) have also presented the effect of initial cyanide 
concentration on gold extraction in the presence of activated carbon. This effect is 
given in Table 5. 



Table 5 The effect of 12 g of activated carbon on gold extraction after 22 h 

of leaching (Rees and van Deventer, 2000). 



Initial free 


Extraction without 


Extraction with 


CN" 


carbon 


carbon 


(ppm) 


(%) 


(%) 


310 


75.0 


84.8 


135 


72.9 


78.4 


115 


68.2 


74.3 


125 


62.6 


73.0 



It can be seen that the level of cyanide concentration has a slight effect on 
extraction in both the presence and absence of activated carbon. The gold 
extraction is enhanced by adding the activated carbon. 
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Rees and van Deventer (2000) further performed gold extraction for different types 
of ore, as shown in Table 6. They observed that the presence of activated carbon 
resulted in an important improvement in the extraction for all ore types. 



Table 6 The effect of 4 g of activated carbon on gold extraction at an initial 

free cyanide concentration of 200 ppm after 24 h (Rees and van 
Deventer, 2000). 



Ore 


Extraction 


Extraction 


Final 




without carbon 


With carbon 


[ON"] 




(%) 


(%) 


(ppm) 


Oxide ore 


66.9 + 2.5 


80.4 + 4 


40 


Pyrite concentrate 


13.2 + 0.2 


16.0 + 0.2 


0 


Copper concentrate 


1.5 + 0.02 


3.2 + 0.16 


0 



The selectivity of activated carbon for gold over other metal cyanide species 
increases when cyanide concentration is increased (Marsden and House, 1992). 
Rees and van Deventer (2000) presented that the optimum level of initial free 
cyanide concentration in the presence of activated carbon is relying on ore type. 
They found that the improvement of gold extraction in the presence of activated 
carbon is not strongly dependent on cyanide concentration for the sulfide and 
copper ores. 
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5 GOLD CYANIDATION PROCESSES 

The adsorption of gold from aqueous solutions onto activated carbon has 
apparently been known since 1847 (Fleming, 1992). The continuous processes of 
gold adsorption, namely the carbon-in-pulp (CIP) and the carbon-in-leach (CIL) 
processes, have been widely used since 1970s (de Andrade Lima, 2007). Marsden 
and House (1992) have presented in Figure 20 that CIP and CIL processes have 
been the dominant gold recovery methods in use worldwide, accounting for 
approximately 44% of world production. Other methods are about 30% for solid- 
liquid separation and zinc precipitation (CCL & Zn) and 19% for flotation and 
gravity concentration (MISC). These adsorption processes are significantly used 
since they avoid the need for solid-liquid separation stages, for example, 
thickening and filtration, thus can reduce the operating costs (Marsden and House, 
1992). 




Figure 20 World gold production by recovery method (Marsden and House, 
1992). 

5.1 Carbon-in-pulp (CIP) Process 

Metal recovery or concentration by adsorption onto activated carbon by the 
carbon-in-pulp (CIP) has become a successful gold extraction technique after 
milled ore is leached by cyanide solution in aerated alkaline slurries. It comprises 
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of two main processes which are firstly the leaching step followed by the carbon 
adsorption as schematically shown in Figure 21. 




Figure 21 Flow sheet of a modern CIP plant (Fleming, 2004). 
5.1.1 CIP Leaching Section 

The gold dissolution is largely controlled by the factors such as cyanide and 
oxygen concentrations, pH and temperature. The high concentrations of oxygen 
and cyanide improve the gold leaching and thus reduce the residence time and 
increase throughput (Ellis and Senanayake, 2004). 

5.1.1.1 Process Description of CIP leaching section 
The ore is first reduced in size (typically 80% passing 75 um) in crushing and/or 
ball milling circuit. After comminution the pulp is normally dilute and thickening 
might be performed to increase the pulp density to about 50% solid by mass. 
Leaching reagents in the form of cyanide and an oxidant such as air or oxygen are 
added after thickening (Stange, 1999). The creation of a large number of small 
oxygen bubbles and dispersing them in the slurry long and deep enough provides 
adequate oxygen concentration for gold dissolution (Ellis and Senanayake, 2004). 
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In a typical carbon-in-pulp (CIP) process, the milled ore is leached in a series of 
agitated vessels, typically about six vessels, each having a retention time of about 
four hours. The agitated vessels comprise of a conventional mechanical agitator 
including blades and shaft. The vessels preferably have a cover for the solution to 
minimize the transfer of oxygen from the slurry to air. The cover can be either a 
stationary cover or a floating cover, such as a disc-shaped cover which has a 
generally flat top surface, and generally concave bottom surface which is actually 
in contact with the slurry. In the agitated leaching reactor the gold is dissolved 
from the pulp. Pure oxygen or natural air, which can be bubbled into the vessels, is 
injected at the lower part of the leach tank (Brison et al., 1989). Recently, de 
Andrade Lima and Hodouin (2006) have investigated the configuration for a gold 
leaching circuit based on gold price and cyanide cost with five equal-sized 
reactors. They concluded that for high gold content, 0.5-1.2 mg/kg, and low 
cyanide concentration, 250-500 mg/dm , the conventional in- series configuration is 
the most recommended. However, for high cyanide concentration, 750-1000 

3 3 

mg/dm , and low gold content less than 0.5 mg/dm , the in-parallel configuration is 
economically the most recommended. 

5.1.2 CIP Adsorption Section 

Activated carbon or activated charcoal is the most widely used as adsorbent due to 
its efficient properties, for example, large specific surface area of approximately 
1000 m /g and strong adsorptive capacity (Marsden and House, 1992; Rees et al, 
2000). Various types of activated carbon can be used. However, it is often 
preferred to use activated granular coconut shell carbon because it is harder than 
most other available types. It is more resistant to breakage and abrasion and has 
large loading capacity (Kunter and Turney, 1980; Coetzee and Grey, 1999; Yalcin 
and Aral, 2002). It is widely accepted that gold adsorption onto activated carbon is 
a complex phenomenon. Many studies have proposed a number of significant 
mechanisms which can be divided into four categories (Marsden and House, 1992; 
Vegter, 1992): 

(a) Adsorption as Au(I) cyanide ion 
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(b) Adsorption as molecular AuCN 

(c) Reduction and adsorption as metallic gold 

(d) Adsorption in association with a metal cation such as Ca 2+ 

The mechanism of gold recovery from cyanide liquors is presented in Equation 
(34). The dissolved gold cyanide is adsorbed as the cation gold cyanide ion pair 
(McDougall, 1985; Marsden and House, 1992; Vegter, 1992): 

M" + + nAu(CN) 2 ~ - M" + [Au(CN) 2 ~] n (34) 

where the ion pair, M" + [Au(CN) 2 ] «, is the adsorbed gold species. M can be H, 
Na, K and Ca. The rate steps are the mass transfer of gold cyanide and cations 
across the solution film surrounding the carbon particles to the outer surface of the 
carbon particles, transfer in the solution filling the pores, the adsorption reaction 
and possibly transfer of the adsorbed cation gold cyanide ion pair along the pore 
walls by surface diffusion (Vegter, 1992). It has been claimed that film transfer 
controls the adsorption rate of aurocyanide onto the activated carbon in the initial 
stages of gold recovery (Coetzee and Gray, 1999). 

The equilibrium of the adsorption reaction is influenced by many factors, for 
example, temperature, gold cyanide concentration, presence of organic solvent, 
type of cations in the solution and the ionic strength of the solution (Vegter, 1992). 
In the initial rate of gold adsorption, gold is adsorbed on the surface of macropores 
and possibly mesopores in activated carbon at a high rate. After that the adsorption 
rate decreases as equilibrium is approached, as shown in Figure 22. 
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Figure 22 An example of kinetics of gold loading onto activated carbon 
(Marsden and House, 1992). 

Marsden and House (1992) have presented a first-order rate equation for gold 
adsorption onto carbon as follows: 

log CauJX) = k ads t + log CZ, L (35) 

where Ca u ,l(^> g°ld concentration in solution at time t, mg/dm 3 
C'2 L initial gold concentration, mg/dm 
kads rate constant for gold adsorption, s" 1 

Ibragimova et al. (2007) have also presented a first-order reaction where the rate 
constant k ads of [Au(CN) 2 ]~ ion adsorption on activated carbon is presented as 
shown in the following equation: 

dCAu.c I dt = kads Cau.l (36) 

where Ca u ,c and Ca u ,l are the gold concentrations in the carbon and solution phases, 
respectively and the adsorption rate constant, k ads , is dependent on a multitude of 
parameters. Coetzee and Grey (1999) have presented a typical adsorption profile 
for the recovery of aurocyanide by the activated carbon in a batch reactor as shown 
in Figure 23, indicating the regions and points of importance. The diminishing 
aurocyanide decay period replaces the constant rate adsorption when the rate of 
intraparticle diffusion becomes the rate controlling factor (Coetzee and Grey, 
1999). 
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Figure 23 A typical adsorption profile of aurocyanide onto activated carbon 
(Coetzee and Grey, 1999). 



The equilibrium isotherm is further illustrated in Figure 24. It is an empirical 
equation relating the concentration of the chemical species in the liquid phase to 
their concentration on the surface of the activated carbon. 
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001 0 10 10 laOcom 
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Figure 24 Equilibrium adsorption isotherm for loading of gold on carbon 
(Marsden and House, 1992). 

An empirically developed equilibrium gold loading capacity (K value) has been 
suggested to be used for the evaluation of carbons used in gold adsorption systems. 
It is obtained by contacting various weights of carbon with a standard borate 
buffer-gold solution for a fixed time. The K value is interpolated as the carbon 
loading in equilibrium with a residual gold solution concentration of 1 mg/dm 
(Marsden and House, 1992). Figure 24 shows the Freundlich isotherm which has 
been used to depict gold adsorption on activated carbon. Ca u ,c represents the gold 
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concentration in carbons, Ca u ,l is the gold concentration in solution, b is the 
Freundlich adsorption constant and n' is the Freundlich exponent. In the case of n' 
= 1, the Freundlich isotherm reduces to linear isotherm which is accurate for low 
concentration systems, as is usually true for this case. The Langmuir isotherm is 
also used to describe gold adsorption (de Andrade Lima, 2007). 

5.1.2.1 Process Description of CIP adsorption section 

After leaching, the pulp is then introduced to the CIP adsorption process. The CIP 
consists of three essential stages (Marsden and House, 1992; de Andrade Lima, 
2007): 

(a) Adsorption: the dissolved gold in the pulp is loaded onto activated carbon 

(b) Elution: gold is removed from carbon into an alkaline cyanide solution 

(c) Electro winning: gold is removed by an electrical process from the alkaline 
cyanide solution and deposited on electrodes. 

The CIP adsorption process occurs continuously in a cascade of large agitated 
tanks, typically approximately six vessels each having a retention time of about 
one hour. Conventional components of the vessel include the mechanical agitator 
including blades and shaft, the slurry inlet and outlet covered by a screen and the 
carbon inlet and outlet. Non-conventional components of the tank also include a 
sparger located adjacent the tank bottom for sparging oxygen into the tank in order 
to enhance adsorption efficiency onto the activated carbon (Brison et al, 1989; 
Marsden and House, 1992). Another non-conventional vessel comprises of top 
sealing in the tank in order to maintain an oxygen atmosphere of about one 
atmosphere pressure, or significantly greater (Brison et al, 1989). The pulp is 
agitated and flows by gravity from one tank to the next by overflow. The solution 
first comes into contact with the activated carbon in the first tank. The fresh or 
reactivated carbon is introduced at the tail end of the process and is 
countercurrently transferred either continuously or periodically with respect to the 
pulp. The carbon in each stage becomes loaded to pseudo-equilibrium, depending 
on the gold concentration in the solution in each stage. In the first stage, the carbon 
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has the highest gold loading and it is contacted with the highest concentration of 
the solution and vice versa. The slurry is discharged through the outlet of the tank 
to a disposal site for the ore tailings. The gold containing carbon is treated by an 
elution step to desorb the metals which are further recovered in electro-winning 
process. The eluted carbon is regenerated with acid to remove the contaminants. 
After that it is recirculated into the adsorption-elution circuit. 




Figure 25 Schematic diagram of a carbon- in-pulp or carbon-in- leach process 
with three tanks, showing inter-stage screens (IS), the screen at the 
exit of the first tank (S), and the carbon transfer pumps (P) (de 
Andrade Lima, 2007). 



The process involves the use of inter-stage screens (IS) where carbon is retained in 
each adsorption stage due to the fact that it is much larger than the ore particles. 
Figure 25 depicts the inter- stage screens at the reactor exit that confine the carbon 
in the tanks until the transfer of pulp is activated. The screens have apertures which 
are slightly smaller than the activated carbon but large enough to allow free flow of 
the slurry between the stages, typically 0.6-0.8 mm (Stange, 1999). The efficiency 
of inter- stage screening has received much attention since it is claimed to be one of 
the biggest challenges in the application of carbon adsorption. 

The activated carbon can be used in a variety of size ranges, usually 1.7 x 3.4 mm, 
1.2 x 3.4 mm and 1.2 x 2.4 mm. The tanks normally contain of 10-25 grams of 
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carbon per liter of pulp (0.5 to 1.2% by volume carbon). Loaded carbon values on 
operating plants range from 300 to 20,000 grams of gold per ton of carbon 
(Marsden and House, 1992). 

5.2 Carbon-in-leach (CIL) Process 

Conventional carbon-in-leach (CIL) process is similar to the CIP process except 
that the leaching and adsorption steps take place in the same tanks simultaneously. 
Carbon is added in the leach tank so that gold can be adsorbed onto carbon almost 
as soon as it is dissolved by cyanide solution. The CIL process is mainly used to 
treat ores containing organic matter and preg-robbing species that are able to 
adsorb [Au(CN) 2 ]~ anions from solution by carbonaceous species in ore, clay-type 
minerals and other mineral species such as graphite (Marsden and House, 1992; 
Rees and van Deventer, 2000; van De venter et al, 2004; de Andrade Lima, 2007; 
Ibragimova et al., 2007). The process has the advantage of lower capital 
investment cost than separate leaching and carbon adsorption systems and it also 
improves gold extraction by the presence of carbon. The CIL process is a good 
alternative, but there are also some disadvantages. For example, it requires a larger 
carbon inventory resulting in a larger in process tie-up of gold. Carbon and gold 
losses due to carbon attrition are higher and carbon loading will be lower than in 
CIP process, meaning that more carbon has to be stripped (Marsden and House, 
1992). 

In a typical CIL technique, the ground and thickened ore slurry typically passes to 
a series of about six agitated vessels, each having a retention time approximately 
four hours (Brison et al, 1989). Figure 26 schematically illustrates the CIL 
process. The activated carbon is mixed with the ore pulp in an agitated leach tank. 
A counter-current flow is maintained between the gold ore slurry, containing gold 
cyanide and oxygen in solution, and the activated carbon. The leached slurry flows 
down the cascade while carbon is transferred periodically upstream. In the first 
tank of the series, leaching of the fresh pulp is a primary activity. In the latter 
tanks, adsorption onto the activated carbon is dominant as the concentration of 
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gold in the solution increases and the fresh carbon is loaded into the system. The 
pulp then passes through the carbon inter-stage screens the same manner as in the 
CIP process shown in Figure 25. The aperture of the inter-stage screens allows the 
leached pulp to pass through by gravity. The loaded carbon then passes to a loaded 
carbon screen where the slurry is screened. The carbon is further stripped and 
thermally regenerated to restore the active surface and to increase adsorption 
activity before reuse. The gold solution moves to electrowinning process where 
gold and other metals are precipitated on steel wool cathodes (Brison et al., 1989; 
U.S. Protection Agency, 1993; van Deventer et al., 2004). 
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Figure 26 The carbon- in- leach process (Outokumpu, 2004). 



According to Brison et al. (1989), typical conditions of the ore slurry would be 
50% solids (-100 mesh), solid specific gravity of 2.7 and slurry specific gravity of 
1.46. They also experimented with oxygen atmosphere, pH 9-11, 0.3 grams of 
NaCN per liter. The results show the gold leaching in the range of 91-92% in six 
hours and low cyanide consumption. 
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5.3 Cyanidation at High Pressure and/or Elevated Temperature 

The direct oxidative pressure cyanidation has been recently presented as a new 
technology where high pressure is a main operating condition. This improved 
technique is applied in autoclaves where oxidation and dissolution of gold can take 
place simultaneously. Several metallurgical plants under projects, for example, 
Goldstrike, McLaughin, Campbell Red Lake, have successfully used pressure 
leaching as pretreatment of refractory gold ores (Fleming, 1992). This high 
pressure leaching can treat pyrite refractory gold concentrations as well. The 
decomposition products of the oxidation of refractory gold ore are dissolved 
species and can be precipitated under certain condition. The partial decomposition 
reactions of refractory gold ore are presented as follows (Fleming, 1992; Parga et 
ah, 2007): 

2FeS 2 + 70 2 + 2H 2 0 -> 2Fe 2+ + 4S0 4 2 " + 4H + (37) 

FeS 2 + l/20 2 + 2H + -> 2Fe 3+ + H 2 0 (38) 

Fe 7 S 8 + 32H 2 0 + 62Fe 3+ -> 69Fe 2+ + 64H + + 8S0 4 2 " (39) 

2FeAsS + 70 2 + 2H + + S0 4 2 " + 2H 2 0 -> 2Fe 3+ + S0 4 2 " + 2H 3 As0 4 (40) 

The reactions taking place at high pressure also provide the formation of elemental 
sulfur layer and thus allowing cyanide and dissolved oxygen to enter the locked 
gold (Parga et al, 2007). 

Recently Lorenzen and Kleingeld (2000) have presented a novel jet reactor 
technology, with specific regard to gas-liquid-solid cyanidation system using gas- 
liquid jet reactor and solid-liquid jet reactor. They also presented gas-liquid jet 
reactor designs which were determined to be the most efficient configuration as 
shown in Figure 27. The design is for bringing gas into contact with liquid, using 
high levels of turbulence and impingement zone. 
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Figure 27 Simplified diagrams of the gas-liquid reactor designs: (A) Kite- 
shaped reactor and (B) c-shaped reactor (Lorenzen and Kleingeld, 
2000). 

The result shows that a gold recovery of 99% was achieved after only 2¥i hours of 
contact time when using solid-liquid reactor under condition of 1 kg KCN/t, pH 10 
and high pressure at 20-30 MPa. They concluded that these new jet reactors 
connected in series can be used in cyanidation process to increase gold recovery 
(Lorenzen and Kleingeld, 2000). 

5.3.1 Operating Conditions 

It has been found that the oxidative pressure cyanidation can be carried out at 
temperature below about 120 °C and at pressure about 1000 kPa (Corrans et al., 
1993; Hourn et al., 1999). Sundkvist (2002) has presented that the high yield of 
gold recovery, 95% approximately, was obtained in several trials at a temperature 
of 60-80 °C. The best condition for gold production has been found at a 
temperature of 75 °C, at pH 11-11.5. Parga et al. (2007) also experimented at 
operating conditions of 60-200°C, 1% NaCN, PbO 100 g/t, 0.6 MPa oxygen 
pressure and at pH 11.2. The results show that increasing the temperature above 
80°C decreases the gold extraction since the oxidation of cyanide ions is too rapid. 
This results in the formation of CO2 gas and ammonium ions. Consequently, at 
these conditions with the residence time of 60 min the gold recovery exceeds 96%. 
Figure 28 illustrates the percentage of gold extraction at operating conditions 
compared to ambient conditions. 
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Figure 28 A comparison of gold extraction at ambient conditions and high 
pressure. Condition: % solid = 20; pH = 11; speed = 300 min" 1 
(Parga et al, 2007). 



5.4 Intensive Cyanidation 

Intensive cyanidation is a proven technology operating at high gold and reagent 
concentrations and/or high temperature and pressure to increase gold dissolution 
rate. The intensive cyanidation is applied for gravity concentrates containing 
coarse gold. Oxygen may be introduced as air, pure oxygen or mixture of the two 
to gain elevated oxygen partial pressure (Marsden and House, 1992). Table 7 
presents the different leaching conditions between the intensive cyanidation and 
conventional one. 



Table 7 Comparison between leaching conditions in intensive and 

conventional cyanidation (Longley et al, 2003). 



Conditions 


Intensive 


Conventional 


Factor 


cyanidation 


cyanidation 


Au, g/t 


20,000 


2 


10,000 


CN, ppm 


20,000 


500 


40 


0 2 , ppm 


15 


8 


2 


Au dissolution rate, g/h/t 


1000 


0.1 


10,000 



The conventional methods such as high-speed agitator in a tank and vat leaching 
have shown problems due either to excessive wear and high energy requirement 
for the agitated systems and loss of fine gold in the vat type systems resulting in 
poor recoveries. Gold recovery on the production tables is relatively low, often 
between 30 % and 60 % (Longley et al., 2003). Many researchers, therefore, have 
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attempted to develop new machines that can very efficiently and economically 
concentrate gold by gravity. Both the ConSep ACACIA Reactor and Gekko 
System's InLine Leach Reactor are widely used in intensive cyanidation, with gold 
recoveries more than 95 % (Mineral Engineering International, 2005). These new 
devices are able to bring down the operating costs as well. 

5.4.1 ACACIA Reactor 

The ACACIA Reactor is widely known as The ConSep ACACIA Reactor, shown 
in Figure 29, where ConSep is the supplier. This reactor is a high efficiency 
leaching device using a high intensity cyanidation method to achieve very high 
recovery. The process consists of an upflow fluidized reactor, which has been 
developed to produce ideal solid-liquid interaction in order to maximize gold 
leaching reaction kinetics without mechanical agitation. The product from the 
process is in the form of cathode gold which is ready for traditional smelting step 
(Campbell and Watson, 2003). 




Figure 29 



The ConSep ACACIA Reactor (Consep Pty Ltd). 



5.4.1.1 Process Description 
Figure 30 presents the basic flow sheet of the ACACIA Reactor. 
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Figure 30 Basic flow sheet of the ACACIA Reactor (Watson and Steward, 
2002). 

The ACACIA is operated in batchwise. The gold concentrate is pre-washed to 
remove suspended solids that maybe entrained in the gravity concentrate and may 
interfere with the subsequent gold electrowinning. The concentrate solids are then 
fluidized for 30 minutes, approximately. The fine solids are introduced into the 
upflow of the leach solution which is pumped at a rate sufficient to fluidize the 
concentrate. The leach solution is mixed in the Reaction Vessel Feed Tank. The 
concentrated cyanide (NaCN), caustic soda and portable water are added into the 
tank. Typically high leach recoveries can be achieved in less than 8 hours. The 
leach cycle must be completed within 16 hours. The pregnant solution is then 
transferred to the main eluent tanks for eletro winning (Watson and Steward, 2002). 

According to the experiment of Watson and Steward (2002) at temperature 
between 50 to 65 °C with H2O2 and non oxygen based oxidant, the results show 
that the gold in the gravity concentrate is recovered over 98% within 8 hour leach 
time at any condition as shown in Figure 31. The forth case provides the highest 
gold recovery, 99.8% under condition of 80 liters of solution, 1600 grams of 
cyanide, 280 grams of NaOH and concentrate grade of 1.5% Au. 
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Figure 31 Recoveries during commission at different conditions in ACACIA 
Reactor (Watson and Steward, 2002). 



5.4.2 InLine Leach Reactor 

The InLine Leach Reactor (ILR) is a commercial, either continuously or batch wise 
operate intensive cyanidation reactor. The batch models of the InLine Leach 
Reactor have been developed to treat higher- grade concentrates between 1000 and 
20,000 g Au/t, while the continuous InLine Leach Reactors have been developed to 
treat lower-grade concentrates above 50 g Au/t. This allows flexible application of 
the intensive cyanide leaching technology for the treatment of different primary 
concentrates. Oxygen gas or peroxide can be used as the oxygen source. Typically, 
the level of dissolved oxygen is above 20 ppm, pH at 13.5 and the cyanide level is 
2% or 20,000 ppm (Katsikaros N., Gekkos System Pty Ltd, private conversation, 
2008). The leach times range from 6 hours up to 24 hours (Gray, 2000). Figure 32 
illustrates the area of application for each type of ILR in relation to concentrate 
mass and Au grade. 
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Figure 32 Areas of application of the continuous and batch ILR (Gekko 
Systems PtyLtd). 

The commercial name for the reactor is Gekko 's System InLine Leach Reactor 
where Gekko is a manufacturer. The ILR works on the principle of laboratory 
bottle roll to keep the solids in contact with the reagents. It consists of a horizontal 
drum rotating at approximately 1 min" 1 . This allows the reactor to agitate and leach 
coarse material, typically up to 6 mm, at high percent solids. The solids are not 
suspended in the reactor as in typical slurry but they are well mixed so that all 
solids and solution are contacted and the reagents are intimately mixed (Katsikaros 
N., Gekkos System Pty Ltd, private conversation, 2008). Figure 33 presents a 
typical of the InLine Leach Reactor. 




Figure 33 InLine Leach Reactor, 



model ILR 2000 (Gray, 2000). 
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5.4.2. 1 Process Description of the Batch InLine Leach Reactor 

The concentrates from primary recovery device are fed to the feed cone for de- 
watering. The solids are stored in the feed cone until the beginning of each leach 
cycle. The reagents are added and stored in the solution storage tank. During the 
leaching step the solution containing the reagents is continually recirculated from 
the solution storage tank to the reactor feed. The solution is then overflown to the 
sump before being returned to the solution storage tank. Figure 34 shows the 
flowsheet of the Batch InLine Leach Reactor. 




Figure 34 Batch InLine Leach Reactor flowsheet (Longley et ai, 2003). 



The reactor drum is horizontally rotated only fast enough to ensure that fresh 
solution is mixed through the solid. The action of the rotating drum draws air into 
the slurry in order to ensure that no portion of the solids can become oxygen 
depleted resulting in no occurrence of gold re-precipitation onto steel particles. 
During drum loading all solids are retained within the leaching system. At the 
completion of leaching the drum is stopped and the solution in the drum is drained 
into the sump and further pumped to the solution storage tank. Some fines particles 
are carried through the drum and recirculated with the solution. The solution is 
then drained from the bottom of the tanks conical base. Wash water is added to the 
drum to wash entrained pregnant solution from the solids and drained out then 
clarified. The leached washed solids are emptied from the drum by reversely 
rotating the drum. The solids then move to the sump while water is continuously 
added to pump back to the mill circuit. The pregnant solution can either be 
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transferred into the dedicated electro winning tank and cell recovery or directly 
added into the elution circuit for recovery in the main plant electrowinning cells 
(Longley et al., 2003). Figure 35 presents the leach and electrowinning 
performance of the Batch InLine Leach Reactor from the gold mine in Indonesia, 
with 10 ppm of oxygen, 24 hours operating cycle. 
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Figure 35 Leach and electrowinning results for Indonesian ILR2000BA 
(Longley et al., 2003). 



The results show very high recoveries with an average of 99.5% gold dissolution in 
the gravity concentrate. Gold electrowinning recoveries averaged 98% on the leach 
solutions, from gold feed grades of approximately 10,000 g/t. Residue grades of 50 
g/t are regularly obtained (Longley et al., 2003). 

5.4.2.2 Process Description of the Continuous InLine Leach Reactor 

The ILR operates continuously in conjunction with an existing or dedicated 
electrowinning cells with the barren solution recycled back to the ILR to minimize 
the use of reagent. The continuous intensive cyanide process including the ILR and 
the dedicated electrowinning cell is shown in Figure 36. 
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Figure 36 Continuous intensive cyanide process (Gray et al, 2003). 



The feed slurry is introduced to the feed cone where it is continuously dewatered. 
The cone holds the solids and water is overflowed to the ILR tails sump. The 
thickened feed is then introduced to the reactor drum. Fresh reagents and return 
barren solution from the eletrowinning circuit are added in the feed. The drum 
rotates only fast enough to ensure that fresh solution is mixed through the solids. 
Solution and solids then move to the pregnant solution sump which is divided into 
three interconnected sections. The first section allows the solids and solution to be 
pumped to the tails settling cone. The solids settle at the cone bottom and are 
thickened then moved to the inclined dewatering screen. The solids are dewatered 
about 83% by weight and dropped from the end of screen into the tailings 
discharge sump. Any underflow solids from this screen are returned to the pregnant 
solution sump. The overflow from the tailings settling cone goes to the second 
section of pregnant section sump. The third section incorporates a set of lamella 
plates to complete the clarification of liquor. The overflow from the lamella plates 
is then pumped to the solid settling vessel to ensure that the minimal solids reach 
the electro winning cell. Flocculent may be used to improve the clarification 
performance. The barren solution overflows from the tail of electrowinning cell 
and is returned back to the feed of reactor drum. The system provides excess 
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solution, which is released as overflow from a small tank located at the reactor feed 
(Gekko System Pty Ltd.). 



Longley (2004) has investigated the gold leach performance using technology 
combined gravity/flotation concentrate prior to ILR intensive cyanide, with 0.5 % 
cyanide concentration and 40 ppm dissolved oxygen levels. This combining 
method was found to be efficient since it decreased capital costs and improved 
safety due to reduced cyanide consumption. The results show that 92.9 % gold 
recovery was obtained after 8 hours and 95.9 % after 24 hours. 
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Figure 37 ILR intensive cyanide leach recovery curve (Longley, 2004). 
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6 GOLD CYANIDATION MODEL 

In this work, a numerical model for gold leaching is developed to improve the 
reliability of the model predictions. The model takes into account the effects of 
surface reaction, internal diffusion and decreasing of active surface area during the 
leaching. A numerical shrinking-core model (SCM) is considered since it is 
flexible to use with multiple reactions and its kinetic model can be developed with 
different reaction orders. Furthermore external mass transfer limitations can be 
easily involved in numerical models. For shrinking-core model, the zone of 
reaction moves into the solid, leaving behind the completely converted material. 
Thus, at any time there exists an unreacted core of material which shrinks in size 
during reaction (Levenspiel, 1999). 

6.1 Assumptions and Approximations for Gold Leaching Model 

Leaching is a heterogeneous reaction which includes the decreasing of surface area 
as the reaction proceeds. In order to develop the model of gold cyanidation, the 
followings simplifying assumptions are taken as below. 

1. The system of gold cyanidation is assumed to be well mixed so that the solid- 
liquid mass transfer resistances can be neglected. 

2. The reactions taking place follow Equation (7) where no side reactions are taken 
into account. 

3. The leaching reaction is assumed to take place inside the spherical particle so the 
shrinking-core model is considered as shown in Figure 38. 

4. The cyanide concentration is in excess when compared to oxygen concentration. 
Therefore, the concentration of cyanide at the reacting surface does not have to be 
solved. 

5. Pseudo steady-state conditions are assumed inside the particle and the particle 
concentration profiles are independent on time. 
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Figure 38 Illustration of the reacting surface and passive layer in the 
shrinking-core model: dp is the thickness of the passive layer. 



6.2 Kinetic Model for Gold Leaching 

In this current study, two leaching rate equations are used to describe the gold 
cyanidation kinetics. The leaching is normally presented as the amount of metal 
extraction over time (the metal extraction curve). The rate equation can be 
expressed with a term for each reactant as follows: 



r _ h.fOC c p C Y 
" ~ CN~ 0 2 ^-Au,s 



2 1 

mol m" s" 



(41) 



where C Au s surface concentration of gold in particle, mol/m 



C cyanide concentration, mol/nr 



C 0 oxygen concentration, mol/m 3 



k overall rate constant, unit depends on reaction orders 



reaction rate at surface, mol/m s 



The reaction of gold is considered as a surface reaction. Gold dissolution begins 
when the gold in the particles reacts with the reagents at the particle surface. The 
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dissolved gold diffuses out of the particle due to solution concentration gradient. 
According to the shrinking-core model as illustrated in Figure 38, a passive layer is 
formed on the particle surface above the reactive surface. The thickness of the 
passive layer will continually increase as the reaction proceeds. Oxygen diffuses 
from solution through this layer to the reactive surface. In the rate equation it may 
be assumed that the gold concentration at the particle surface is constant at any 
time since it is characteristic to the concentrate. Consequently, the reaction rate 
constant, k, and gold concentration, Ca u ,s can be lumped, as given: 

r s =KC a CN .C^ s molmV (42) 

"3 

where C 0 ^ s oxygen concentration at the reacting surface, mol/nr 

k\ overall reaction rate constant, unit depends on reaction orders 

Another kinetic equation is obtained if the reaction rate is assumed to depend also 
on the gold concentration in the particles. Some amount of gold in the particles 
remains unreacted thus the non-reacting gold should be removed from the total 
gold concentration in the model. The gold concentration is presented as the 
difference of gold concentrations between its concentration and its residual gold 
concentration after infinite leaching time. The rate equation becomes: 

r s = k 2 C a CN _C^ s (C Au -CZY mol in 2 s" 1 (43) 

where C Au gold concentration in particles, mol/m R 3 

C°° ku gold concentration in particles after infinite leaching time, mol/m R 3 
k 2 overall reaction rate constant, unit depends on reaction orders 

The gold dissolution rate can be given as: 

mol m' 3 s" 1 (44) 
where A, surface area of one particle during leaching, m 2 



dC 



Au 



dt 



= - r sA n 
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"3 

n number of particles per volume, l/m R 



The surface area of sphere particle related to shrinking-core model can be 
calculated as follow: 



A i =71 d p 



( r \% 

^ Au 

y Aujini j 



(45) 



where Ca u ,m initial gold concentration in particles, mol/m R 3 
dp particle diameter, m 

6.3 Internal Mass Transfer of Oxygen in Gold Leaching 

Mass diffusion rates within the porous residue are difficult to assess. This is 
because it is impossible to define the shape of the channels through which the 
transfer takes place. However, the rate of the transfer in a spherical diffusion layer 
can be approximately obtained. 



C 




Figure 39 Illustration of a spherical diffusion layer. 

The oxygen concentration at the reacting surface can be solved numerically from 
mass balance under assumption of pseudo-steady state inside the particle. The 
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following equations present the concentration profile in a spherical diffusion layer 
as illustrated in Figure 39. 

Cr(r o - r t ) = C t r. (r o - r) + C o r o (#■-/;.) (46) 



C(r) = CMr 0 -r) + C 0 r 0 (r- ri ) (4?) 
r(r a -r t ) 

where C concentration in spherical diffusion layer, mol/m L 
d concentration at inner passive layer, mol/m L 3 

-i 

C„ concentration at outer passive layer, mol/m L 

r radial coordinate, m 

r, inner radius of passive layer, m 

r 0 outer radius of passive layer, m 



The concentration gradient in the layer is then obtained as follows: 

dC _ r i r 0 (C 0 -C i ) 
dr r 2 (r 0 -r,) 



(48) 



The equation for mass diffusion rate of oxygen at the inner and outer surfaces can 
be given in terms of concentration gradient and effective diffusivity (D e ) as shown 
below: 



I S~l 

No, =-D e A r -— mols 1 (49) 

dr 

= _ D A r ° (C °i ~ C °^ mo i s -i ( 50) 

' r Mo- r i) 

= _ D A r ^ C o 2 ~ C o 2 ,) ^ s -l (51) 

" r o( r o- r i) 

where A r reacting surface area of particles, m 
A r inner surface area of passive layer, m 2 
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2 

outer surface area of passive layer, m 



2 

D e effective diffusivity in particles, m /s 



No 2 diffusion rate of oxygen, mol/s 

The radius of the reactive surface r, is continuously recalculated as the reaction 
proceeds and gold is leached. According to the shrinking-core model, the reactive 
surface decreases as the reactive layer penetrates deeper and deeper into the 
concentrate particle. Therefore, the active area in Equation (44) continually 
reduces as the reaction proceeds. 

The steady state diffusing rate is equal to the rate of reaction of oxygen within the 
reaction zone. Therefore, the mass balance for oxygen at the reacting surface can 
be obtained and oxygen concentration at ore surface can be solved for Equations 
(42) and (43) as presented, respectively: 



0.25k, c a CN _ c^4-=£>A 



r i(Cp 2 ~ C o 2 , s ) 

r o( r o- r i) 



(52) 



and 0.25k 2 C a CN _ C^ s (C Au -C~ Au Y A,. = D e A, o 



ri(C 02 -C 0iS ) 
r o( r o- r i) 



(53) 



where the value of 0.25 is the stoichiometric constant of oxygen referring to 
Equation (7). 
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7 MODEL PARAMETER ESTIMATION 

The gold dissolution equations which were used to fit the model parameters to the 
experimental data are presented as follows: 

^ = -^C^C^A ( n molmV (54) 

and ^ = -k 2 C^C^ s (C Au -C: u yA in mol m" 3 s" 1 (55) 

7.1 Experimental Data Available 

The articles presented by Ling et al. (1996) and de Andrade Lima and Hodouin 
(2005) are used as sources for the experimental data in order to validate the studied 
models in this work. The experimental variables in Ling et al. (1996) were particle 
size, dissolved oxygen, cyanide and gold concentrations, as shown in Table 8. The 
temperature and pH were assumed to be constant at 20°C and 1 1 , respectively. The 
experiments were carried out in a 1 liter batch-type glass beaker reactor equipped 
with four Teflon coated baffles and a flat-blade turbine agitator. The agitation 
speed was maintained at 750 min" 1 . The amount of ore was 600 g and the total 
slurry volume was around 800 ml. The dominant compositions of the ore are Na- 
feldspar (NaAlSi 3 0 8 ), calcite (CaC0 3 ) and quartz (Si0 2 ). 



Table 8 Experimental conditions in the data sets from Ling et al. (1996). 



Data No. 


Mean Diameter, um 


C CN 'PP m 


C 0i , ppm 




1 


16.51 


104 


8.5 


4.70 


2 


37.42 


104 


8.5 


8.00 


3 


37.01 


13 


8.5 


6.47 


4 


37.01 


25 


8.5 


6.47 


5 


37.01 


104 


8.5 


6.47 


6 


37.01 


156 


8.5 


6.47 


7 


37.01 


104 


20 


6.47 


8 


37.01 


104 


30 


6.47 


9 


46.47 


104 


8.5 


13.20 
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Table 9 Experimental conditions in the data sets from de Andrade Lima and 

Hodouin (2005). 



Data 


Mean Particle Diameter, 

1 1 in 

JLIII1 


C cn 'PP m 


C 0i , ppm 


C Au , ppm 


1 


30 


260 


8 


2.2 


2 


30 


650 


40 


2.2 


3 


30 


260 


8 


2.2 


4 


30 


650 


40 


2.2 


5 


44 


260 


8 


1.9 


6 


44 


650 


40 


1.9 


7 


44 


260 


8 


1.9 


8 


44 


650 


40 


1.9 


9 


63 


260 


8 


1.5 


10 


63 


650 


40 


1.5 


11 


63 


260 


8 


1.5 


12 


63 


650 


40 


1.5 


13 


88 


260 


8 


2.3 


14 


88 


650 


40 


2.3 


15 


88 


260 


8 


2.3 


16 


88 


650 


40 


2.3 


17 


125 


260 


8 


2.2 


18 


125 


650 


40 


2.2 


19 


125 


260 


8 


2.2 


20 


125 


650 


40 


2.2 


21 


177 


260 


8 


1.7 


22 


177 


650 


40 


1.7 


23 


177 


260 


8 


1.7 


24 


177 


650 


40 


1.7 



Table 9 presents the conditions of the experimental data from de Andrade Lima 
and Hodouin (2005). The experimental variables were the dissolved oxygen 
concentration, cyanide concentration and particle mean size as well. The 
temperature and pH were set to be constant at 20°C and 12, respectively. The 
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experiments were carried out in a closed batch reactor of 1000 mL plastic beaker 
equipped with four baffles and four-blade turbine, in which 500 g of ore was 
added. The dominant components of the ore are quartz (Si02), clinochlore 
((Mg,Fe,Al) 6 (Si,Al) 4 Oio(OH) 8 and illite ((K,H 3 0)Al 2 Si 3 A10io(OH) 2 ). 

7.2 Results and Discussions 

The kinetic equation parameters were estimated by the Modest computer software 
(Haario, 2002) using the least squares method. The squared difference between the 
simulations and the experiments were minimized. The estimated parameters were 
obtained by fitting the model to the data. For the model presented in Equation (54), 
the estimated parameters consist of ku a, /? and D e , while parameters &2 and y are 
also estimated in the model of Equation (55). The estimated parameter values 
obtained from Equation (54) using experimental data shown in Table 8 are given 
below in Table 10. 



Table 10 The estimated parameter values and their standard error from 
Equation (54) 



Data Source 


Parameter 


Values 


Std. Error, % 


Ling et al. (1996) 


ki (m 3 /mol) 235 h' 1 


0.555 x 10" 5 


2603.8 


Explained: 9.68% 


a 


2.29 


27.1 




ft 


1.06 


1838.9 




D e (m 2 /h) 


0.101 x 10" 8 


1147.2 


de Andrade Lima 


MmoVm 3 ) 0 - 668 !!- 1 


0.1 14x 10" 6 


39.7 


and Hodouin (2005) 


a 


0.32 


44.4 


Explained: 1.14% 


ft 


0.12x 10" 1 


703.6 




D e (m 2 /h } 


0.610x 10" 7 


90475.6 



The rate equation for gold cyanidation becomes: 

r s = 0.555X10" 5 C^ 9 C^ s mol m 2 h" 1 (56) 
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and r s = 0. 1 14 x 10" 6 C™1 C°; 01 s 2 mol m" 2 h" 1 (57) 

Figure 40 demonstrates the fitting of the gold leaching model to the experimental 
data in both cases. The experimental data are given in Appendix IA. The fitted 
curves present the gold conversion on y-axis at any leaching time (h) on x-axis. 
The gold conversion can be calculated using: 

C -C 

s~\ • Au.ini An ✓ cn\ 

Conversion = (58) 

C 

Au,ini 

where Cau,m 1S m e gold concentration in the solids at initial state (t = 0). 



In Figure 40, the subfigures present the data sets from left to right starting from the 
upper-left corner. The coefficient of explanation of Figure 40 a) is 9.68%, meaning 
that the predicted values are not in good agreement with the experimental data. 
Figure 40 a) illustrates the overshooting taken place in the data sets no. 1, 5, 6, 7 
and 8. The data sets no. 1 and 5 denote the experiments with small particle size 
compared to data sets no. 2 and 9. The cyanide concentration in the data set no. 6 is 
higher than in the other data sets. The data sets no. 7 and 8 have higher oxygen 
concentration than the others. This level of oxygen concentration in the solution 
might be in excess with respect to gold. In Figure 40 b), the coefficient of 
determination is only 1.14%. The prediction curves of the first twelve data sets 
containing small particle sizes indicate that there is some overshooting. These 
curves do not fit well to the experimental data, while the other curves containing 
larger size of particles show a better fit. 
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Figure 40 Gold conversion vs. time (h) for Equation (54); predicted 
conversion compared to experimental data; a) data from Ling et al. 
(1996); b) data from de Andrade Lima and Hodouin (2005). 
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The estimated parameters obtained for the gold leaching model shown in Equation 
(55) consist of k 2 , a, /?, y and D e for both data sources. These estimated parameters 
are given in Table 11. 



Table 11 The estimated parameter values and their standard error for 
Equation (55). 



Data Source 


Parameter 


Values 


Std. Error, % 


Linger al. (1996) 


fc 2 (m 3 /mol) 2 - 73 V 


0.706 x 10" 3 


63.9 


Explained: 81.13% 


a 


1.38 


9.4 




ft 


0.161 


149.0 




y 


2.19 


12.0 




D e (m 2 /h) 


0.1 x 10" 7 


416.8 


de Andrade Lima 


fc 2 (m 3 /mol) 201 V 


0.624 x 10" 1 


167.6 


and Hodouin (2005) 


a 


0.705 


36.2 


Explained: 53.75% 


ft 


0.16 x 10" 1 


1012.4 




y 


2.29 


12.2 




D e (m 2 /h } 


0.515 x 10" 9 


2.7 x 10 5 



Therefore, the rate equation of gold cyanidation becomes: 

r s = 0.706 xl0- 3 C£ C^' (C Au -QJ 219 molmV (59) 

and r s = 0.624 xlO" 1 ^ 05 C°™ (C Au -QJ 229 molm^h" 1 (60) 

Note that the unit for reagents and gold concentrations is mol/m . The comparison 
of the model prediction and the experimental data is presented in Figure 41 for 
both data cases. The experimental data are given in Appendix IB. The coefficients 
of determination for Figure 41 a) and b) are 81.13% and 53.75%, respectively. The 
reaction order for the gold is approximately two, one for cyanide and between 
0.0161-0.16 for the dissolved oxygen. The predicted curves fit to the experimental 
data quite well for both data cases and are also in obviously better agreement with 
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the data, when compared to the previous model, at any concentration of reagents 
and particle size. This implies that the gold leaching also depends on the gold 
concentration in the particles. Therefore, the term of the gold concentration 
difference between its current concentration and the residual concentration should 
not be omitted. The results also show that the gold leaching also depends on both 
cyanide and oxygen concentrations. Therefore, the form of Equation (25), should 
not be considered as the gold leaching rate equation. In addition, the obtained rate 
equation implies that the gold leaching rate is more sensitive to the cyanide 
concentration than oxygen concentration. In practice, it is not possible to 
completely leach all the gold from the particles which means that there is some 
amount of gold remained in the particles. Therefore, the gold concentration term in 
the model should not be considered as a constant. The rate equations, as shown in 
Equations (59) and (60) together with the mass balance equations for oxygen at the 
reacting surface as presented in Equations (52) and (53), show that the continually 
reduced reacting surface area has important effect on gold leaching. Therefore, the 
rate equations assuming the constant surface area, as presented in Equations (22) - 
(24), and rate equation without surface area consideration, as shown in Equation 
(26) should not be used to explain the gold leaching. Furthermore, the rate equation 
obtained in this work is not complicated as presented in Equation (28). 

It is important that the developed models should be as reliable as possible. In a 
reliable model the estimated model parameters should be well identified without 
any severe cross-correlation between the parameters. The reliability of the 
estimated parameter can be analyzed by the Marcov chain Monte Carlo or MCMC 
methods. These methods use the numerical technique in Bayesian statistical 
analysis to analyze the nonlinear model. The MCMC methods present the 
probability distributions of the predicted parameters. 

The reliability of the estimated parameters for both data cases is presented in 
Figure 42. The obtained effective diffusivity coefficient in Figure 42 is multiplied 
by internal scaling factor lxlO" 6 . All parameters in Figure 42 a) and b) are almost 
well identified and there is no any cross-correlation between parameters. The 
posterior probability density distributions in Figure 42 a) appear sharp and the 
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distributions position is about in the center. Parameters k 2 , a and y in Figure 42 a) 
and parameter y in Figure 42 b) have a narrow probability distribution. This can be 
also seen from Table 1 1 where the standard deviation value of these parameters is 
lower than the others. The parameters ft in Figure 42 a) and k 2 , a and /? in Figure 42 
b) are acceptable. The parameter D e is not well- identified and shows relatively 
wide probability distribution increasing the uncertainty and standard deviation and 
thus affect to the accuracy of model. This is due to the fact that internal diffusion 
has no effect in gold leaching. 

The significant of the internal diffusion to the gold leaching is also studied in this 
work. Figure 43 demonstrates the effect of internal diffusion as the ratio of the gold 
leaching rate with internal diffusion over the leaching rate without consideration of 
internal diffusion. Figure 43 a) shows that the curves steeply decrease at the 
beginning of leaching time and slightly increase later. From Figure 43 a), it could 
be explained that gold dissolves rapidly and forms passive layer fast in the 
beginning. After that the thickness of the passive layer is reduced showing the 
effect of internal diffusion becoming less and less important. It might be because 
there is only a small amount of gold on the surface to be leached resulting in 
thinner passive layer until the reaction has no internal diffusion effect (the ratio 
equals to 1). Figure 43 b) shows that the ratio is almost equal to 1 during the whole 
leaching time. This shows that there is no internal diffusion effect to be considered 
in the model. Therefore, the gold leaching is kinetically controlled by surface 
reaction. 
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Figure 41 Gold conversion vs. time (h) for Equation (55); predicted 
conversion compared to experimental data; a) data from Ling et al. 
(1996); b) data from de Andrade Lima and Hodouin (2005). 
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a 




a 





Figure 42 The two-dimensional probability distributions of the estimated 
parameters in Equation (55) from the MCMC analysis; a) data from 
Ling et al. (1996); b) data from de Andrade Lima and Hodouin 
(2005). 
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Figure 43 Effect of internal diffusion in Equation (55). Ratio of the actual gold 
leaching rate with / without internal diffusion consideration; a) data 
from Ling et al. (1996); b) data from de Andrade Lima and Hodouin 
(2005). 
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8 CONCLUSIONS 

In this work the effect of leaching parameters on gold cyanidation has been 
studied. It is found that cyanide and dissolved oxygen concentration have a 
dominant effect on gold leaching. Increasing the reagents concentrations can 
improve gold dissolution rate. Small particle size plays a significant role to gold 
dissolution since it increases the area for surface reaction and the contacting area 
between solute and solution. Based on literature, the presence of other minerals, 
especially sulfide mineral, can retard the gold dissolution since they increase 
cyanide consumption. The oxidative processes can be used to treat them before 
leaching. Other parameters are found to slightly affect the gold cyanidation. 

The parameters of the gold leaching models studied in this work were estimated by 
using the Modest computer software. A numerical shrinking-core model, together 
with simplifying approximations and assumptions, were developed for the 
simulation of gold leaching. The leaching parameters are the dissolved oxygen and 
cyanide concentrations and the difference of the gold concentration in the ore and 
the residual concentration in the ore after infinite leaching. The results show the 
advantage of using gold leaching model based on these parameters since the 
models take into account the effects of the reagent concentrations and the 
unreacted gold while the reduced reacting surface area is involved as well. 
However, the studied models have some limitations since side reactions, 
temperature and solid-liquid mass transfer resistances are not considered. The 
leaching parameters were estimated and fitted to the experimental data obtained 
from two different sources. The reliability of these parameters was studied by 
using the MCMC methods. The reliability can be seen from the cross-correlations 
and distributions of parameter and the standard deviations. The large standard error 
shows lacking sufficient reliability. The parameters a, /?, y and &2 are more reliable 
than D e in the case using data from Ling et al, while y is the most reliable in 
another case. The effect of internal diffusion as the ratio of the gold leaching rate 
with internal diffusion over the leaching rate without consideration of internal 
diffusion was also studied. The result shows that there is no internal diffusion 
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limitation, therefore the reaction is controlled by surface reaction. However, the 
reaction also depends on the gold concentration in the particles, as proposed. The 
reactive surface area is continuously reduced during leaching. Therefore the model 
considering the change of surface area is more accurate than the model assuming 
constant surface area. The rate equations including the estimated parameters are 
shown as follows: 

r s = 0.706 x 10- 3 C£ C™ (C Au - C~ Au ) 2 ' 19 mol m" 2 h" 1 

and r s = 0.624 xlO" 6 ^ 05 (C Au - C" ) 23 molmV 

The rate equations show that the reaction orders are approximately two for gold, 
one for cyanide and a smaller value for oxygen. These reaction orders are at the 
same magnitude as those reported by Ling et al. (1996) and de Andrade Lima and 
Hodouin (2005). From the above equations, it is recognized that a unified kinetic 
model is difficult to establish since the nature of gold ores is complex. However, 
the proposed equations provide reasonable explanation. It can be concluded that 
the rate equation following Equation (43) is more reliable than other models 
presented in literature. This is because the reagents concentration, the particle size 
and the unreacted gold concentration are taken into account. Therefore the 
optimization of gold leaching should further be developed based on this equation. 
Additionally, more experimental data should be obtained in order to improve the 
model fitting. 

For the future works, experiments should be done in stirred tank reactors under 
varying concentrations of reagents at constant temperature and pH in order to 
validate the model obtained in this work. Experiments should be also carried out in 
conditions where the effect of gas-liquid and solid-liquid mass transfer could be 
studied. Various types of ore containing different components of metal ions should 
be studied with this model as well. The latter study can be done under the optimal 
concentration of reagents at high pressure and/or high temperature in an autoclave 
in order to examine the effect of high pressure and temperature to the gold 
leaching. 
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